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TERMINOLOGY 

A  notch  filter's  output  frequency  spectrum  consists  of 
a  single  minimum  in  its  frequency  response.   Usually  it  is 
required  that  the  notch  amplitude  be  zero;  this  leads  to  a 
zero  notch  filter.   Two  further  characterizations  are: 

(a)  An  equal  amplitudes  notch  filter  has  equal  amplitudes 
at  zero  frequency  and  infinite  frequency. 

(b)  An  unequal  amplitudes  notch  filter  has  unequal  ampli- 
tudes at  zero  and  infinite  frequencies  in  its  frequency  re- 
sponse curve. 


INTRODUCTION 

Notch  networks  are  important  in  connection  with  feedback 
amplifier  problems.   In  the  past,  several  notch  networks  with 
equal  amplitudes  and  unequal  amplitudes  at  zero  and  infinite 
frequencies  have  been  considered.   The  purpose  of  this  paper  is 
to  generalize  networks  of  the  Goldman  (I4J  type  and  determine 
conditions  for  existence  of  the  notch  frequencies. 

A  GENERAL  NOTCH  FILTER  OP  THE  GOLDMAN  TYPE 

A  general  notch  filter  resistance-capacitance  network  shown 
in  Pig.  1  is  bisected  into  half-sections,  which  are  shown  in 
Pig.  2  and  Pig.  3. 

If  A  is  the  short-circuit  input  impedance  of  the  half- 
section  and  if  B  is  the  open-circuit  input  impedance  of  the 
half-section,  then  subsequent  calculations  yield 

RCrs2  +  (RC  +  RCr  +  r)  s  +  (r  +  1) 

A  = (1) 

(RCr  +  rCjs^  +  (RC  +  rC  +  C  +  r) s  +  1 


1    s  +  1 

B  =  1  +  -  =  (2) 

s     s 


Bartlett's  representation  theorem  for  symmetric  networks 
yields  the  voltage  transfer  function  T(s).   Actual  calculation 
yields 


T(s)  =  

B  +  A 


[rCs3  +  (2rC  +  C)s2  +  (RC  +  rC  +  C)  s  +  l} 
[(2RCr  +  rC)s3  +  (2RCr  +  2RC  +  2rC  +  2r  +  C)s2 
+  (RC  +  rC  +  C  +  2r  +  2) s  +  l] 


EQUAL  AMPLITUDES  NOTCH  FILTER  OF  THE 
GOLDMAN  TYPE 


Referring  to  Eq.  (3)  of  the  general  notch  filter  with 
R  =  0,  one  obtains  the  transfer  function  T-j_(s). 

Bl  "  Al 
T-,(s)  =  

Bl  +  Al 

rCs3  +  (2rC  +  C)s2  +  (rC  +  C)s  +  1 
rCs^  +  (2rC  +  2r  +  C)s2  +  (rC  +  C  +  2r  +  2)  s  +  1 

Nx(s) 


(3) 


D,(s) 


(W 


Now,  choose  r  and  C  such  that  the  numerator  polynomial  N,(s)  has 
(s2  +  0)2)  as  a  factor.   The  Routh  array  is  applied  to  determine 
values  of  r  and  C.   First  of  all,  consider 

N]_(s)  =  a1s^  +  b-|_s   +  C]_s  +  d]_  (5) 

The  Routh  array  formed  from  the  even  and  odd  polynomials 
of  N1(s)  is 


Jlul 


In  order  to  have  a  common  factor  of  the  form  ( s^  +  or) ,  a  row 
of  zeros  is  required.   The  condition  that  elements  of  the 
third  row  be  zero  is 


;ldl 


=  0 


(6) 


Returning  to  the  original  numerator  polynomial  N^a)  = 

rCs^  +  (2rC  +  C)s2  +  (rC  +  C)s  +  1,  and  identifying  with  Eq.  (5). 

one  obtains 

rC  •  1 
rC  +  C =  0 


2rC  +  C 


C  = 


(7) 


2r^+  3r  +  1 


Equation  (7)  is  one  equation  involving  two  unknown  param- 
eters.  Another  equation,  of  the  form  f(r)  =  0,  is  obtained  by 
requiring  that  the  value  of  C  be  a  maximum.   It  is  clear,  with- 
out loss  of  generality,  that  the  maximum  value  of  C  is  desir- 
able because  of  the  low  sensitivity  of  these  two  parameters  in 
the  neighborhood  of  a  maximum. 

In  order  to  find  the  maximum  C,  one  must  calculate  the 
first  derivative  of  Eq.  (7)  with  respect  to  r. 


dC    2r2  +  3r  +  1  -  r(lj.r  +  3) 

—  =  (8) 

dr       (2r2  +  3r  +  l)2 

Setting  Eq.  (8)  equal  to  zero,  one  obtains 

-  2r2  +  1  =  0 


1 
r  =  -=   ohm  ( 9 ) 

i 

Substitution  of   r  =  j=   into  Eq.  (7)   obtains 
/  2 


C=3-2  ■  J  2   =   0.17158  farad  (10) 

The  notch  frequency  can  be  determined  by 

bxs2  +1=0  (11) 

That  is, 

(2rC  +  C) s2  +  1  =  0 


12 


C(2r  +  1) 


r 

Substitution  of  C  =  in  Eq.  (7)  and 

(r  +  1) (2r  +  1) 

1 
r  =  - —  ohm  in  Eq.  (9)  into  Eq.  (12)  yields 


?   r  +  X 

u2  =  =  2.ij.ll4.2 


u  =  1-5537  rad./sec.  (13) 

Construction  of  a  Routh  array  from  numerator  and  denomi- 
nator polynomials  yields: 

The  Routh  Array 


rC 

2rC    +    2r 

+ 

c 

rC   +   C   +   2r  +   2 

1 

rC 

2rC   +   C 

rC   +   C 

0 

2r 

2r  +   2 

1 

rC 

rC   +   C/2 

0 

1 

1 

C/2 

0 

There  is  no  zero  row,  and  hence  there  exists  no  common  factors 
in  the  numerator  and  denominator  polynomials. 

After  substituting  these  specific  values  of  r  and  C  into 
Eq.  ([).),  the  voltage  transfer  function  becomes 

0.121315  s3  +  0 . i4.H4.2i  s2  +  0.292895  s  +  1 

T-,(s)  =  5 5 

1      0.121315  sJ+  1.82814.2   s2  +  3.707105  s  +  1 

In  order  to  calculate  the  sinusoidal  spectrum,  let  s  =  ju 
and  obtain 


Tj^Ju) 


(1  -  O.I4.H4.2IC02)  +  j  (0.292895"  -  0.121315^) 
(1  -  1.828l4.2u2)  +  j(  3. 707105"  -  0.1213l5u3) 


Prom  this,  one  obtains 

,2    (1  -  0.1+1J+210)2)   +  (0.292895"  -  0.1213l5u3) 
(1  -  1.8281i2a)2)   +  (3.707105"  -  0.121315"3) 


By  inspection,  one  can  see  that 
i       ,2 


I'llJB) 


at  co  =  0 


Tx(  jco)  |=1        at   to  =°^ 

The  log-log  plot  of  T^(jco)    versus  to  is  shown  in  Pig.  7- 
This  frequency  response  shows  that  the  transfer  function  T^(s) 
represents  a  notch  network  with  equal  amplitudes  at  zero  and 
infinite  frequencies.   Data  for  this  plot  are  obtained  with 
Program  2  in  Appendix  C. 


R              l/Cs           l/Cs  H 

vVW\ 1  I        1 I  1 A/VV^ 


r  1 


i/s  -«- 


■J I 


i/a 


Bisection 


Pig.  1.   A  general  notch  filter-  of  the 
Goldman  type. 


A    >■ 


Fig.  2.   Short-circuited  half-section  of  the 
RC  network  in  Pig.  1. 


Vsz: 


Fig.  3-   Open-circuited  half-section  of  the 
RC  network  in  Pig.  1. 


Bisection 


Pig.  Ij..   Goldman  type  resistance-capacitance 
notch  filter. 
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Ai  — *■    ZZVcs 


Pig.  5.   Short-circuited  half-section  of  the 
RC  network  in  Fig.  Ij.. 


l/slt 


Fig.    6.      Open-circuited  half-section    of   the 
RC   network   in  Fig.    I+. 
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UNEQUAL  AMPLITUDES  NOTCH  FILTER 

From  the  General  Notch  Filter,  using  the  same  networks 
shown  in  Figs.  1,  2,  and  3,  one  obtains  the  transfer  function 
as  shown  in  Eq.  (1),  i.e., 


Ac 


To(s   = 


B2  +  A2 

[rCs3  +    (2rC   +   C)s2  +    (RC   +   rC   +   C)s   +   l] 
[(2RCr   +   rC)s3   +    (2RCr  +   2RC    +   2rC   +   2r  +  C)s2 
+(RC+rC+C+2r+2)s+l] 

N2(s) 
D2(s) 


(11+) 


The  Routh  array  of  the  numerator  polynomial's  even  and 
odd  polynomials  is 


RC  +  rC  +  C 
1 


rC 

2rC  +  C 

rC 

RC  +  rC  +  C 0 

2rC  +  C 

In  order  to  have  the  factor  of  (s  +  to  )  ,  the  third  row  should 
be  equal  to  a  zero  row.   Therefore 

rC 


RC  +  rC  +  C  - 


2rC  +  C 


C  = 


(R  +  r  +  1)  (2r  +  I) 


(15) 


13 


The  notch  frequency  co   can  be  obtained  from  the  second 
row  of  the  Routh  array 

(2rC  +  C) s2+  1  =  0 


:i6) 


Vc(2r  +    1) 


Equation  (16)  shows  that  u  is  independent  of  R. 
Rearranging  Eq.  {1%) ,    one  obtains 

C= (17) 

2r2  +  (3  +  2R)r  +  (1  +  R) 

Taking  the  first  derivative  of  C  with  respect  to  r  in  Eq.  (17) 
yields 

dC    f2r2  +  (3  +  2R)r  +  (l  +  R)j  -  r  jV  +  (3  +  2R)! 

—  =  i , '—^ (18) 

dr  |_2r2  +  (3  +  2R)r  +  (1  +  R)j 

For  the  reason  described  previously,  the  maximum  value  of  C 
is  desirable.   This  is  achieved  by  setting  Eq.  (18)  equal  to 
zero,  namely, 


2  v 


r  n 

■    (3  +  2R)r  +  (1  +  R)  -  r[l+r  +  (3  +  2R)J  =  0 
or 

-2r2  +  (1  +  R)  =  0  (19) 

R  =  2r2  -  1  >  0  (20) 

Equation  (20)  implies  that 


11+ 


fT 

r  >  ohm  (21) 

2 


Substituting  Eq.  (20)  into  Eq.  (17)  yields 


2r2  +  (3  +  kr2    -    2)r  +  2r2   l|r2  +  lp?  +  1 


Thus  substitution  of  Eq.  (22)  into  Eq.  (16)  yields 
1 


(22) 


-/c(2r  +  1) 
1 


V2r  +  1  (23) 


2r  +  1 


l+r2+  Ip?  +  1 

First,  assume  co  =  1  rad./sec.   Prom  Eq.  (23)  one  obtains 

1-1 

=  0  ohm                              ( 2i(.) 


2 

Equation  (2l|)  contradicts  the  condition  of  Eq.  (21). 

Second,  assume  u  =  2  rad./sec.   Then  one  obtains  from 
Eq.  (23) 


1.5  ohms. 


2 

This  value  does  satisfy  Eq.  (21). 

Substitution  of  r  =  1.5  ohms  into  Eq.  (20)  and  Eq.  (22) 
yields  respectively 

R=2r-  1=3.5  ohms 


15 


C  =  — =  0.0625  farad 

l\.r   +  lp?  +  1 

Returning  to  Eq.  (11;)  for  specific  values  of  r,  R  and  C, 
one  obtains 

IrCs3  +  (2rC  +  C)s2  +  (RC  +  rC  +  C)  s  +  ll 

T2(s)  =  -rr-± = 

[(2RCr  +  rC)s3  +  (2RCr  +  2RC  +  2rC  +  2r  +  C)s2 

+  (RC  +  rC  +  C  +  2r  +  2)s  +  ij 

1.5  s3  +  i|.0  s2  +  6.0  s  +  16 
12  s3  +  69.5  s2  +  86  s  +  16 

In  order  to  calculate  the  sinusoidal  spectrum,  let  s  =  jco  and 
obtain 

(16  -  I4..0  co2)  +  j(6.0  co  -  1.5  co3) 

T2(jco)  =  

(16  -  69.5  co2)  +  j(86  co  -  12  co3) 

Prom   this,    one   obtains 


|t2(  jco)|2 


(16    -   1:.0   co2)2   +    (6.0   co   -    1.5  co3)2 
(16    -   69.5  co2)2   +    (86   co   -    12   co3)2 

2.25  co6    -    2  oA  -    92   co2   +   256 
lkk  co6   +   2766.25  cA  +   5172   CO2   +   256 


By   inspection,    one    can    see    that 

I  I  2 

T2(  jco)         =   1  at      co  =   0 


JT2(  jco)  |2   =  0.015625  at      co   =  <*=> 
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A  graph  of  [  T2  (  jco)    versus  to  on  log-log  paper  Is  shown 
In  Fig.  8.   This  frequency  response  shows  that  the  transfer 
function  T2(s)  has  a  notch  frequency  and  unequal  amplitudes  at 
zero  and  infinite  frequencies.   Data  for  this  graph  are  ob- 
tained with  Program  3  which  is  given  in  Appendix  C. 
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3 


3 

C 


!q.XOA   ux   q.ndq.no 


AN  EQUAL  AMPLITUDES  NOTCH  FILTER 

A  general  notch  filter  of  equal  amplitudes  and  unequal 
amplitudes  has  been  investigated  in  the  previous  section.   In 
this  section,  the  more  complicated  notch  filter  resistance- 
capacitance  network  shown  in  Pig.  9  is  investigated. 

The  mid  short-circuit  impedance  and  the  mid  open-circuit 
impedances  are  shown  respectively  in  Pig.  10  and  Pig.  11.   Let 
A  be  the  mid  short-circuit  input  impedance;  and  let  B  be  the 
mid  open-circuit  input  impedance.   After  applying  series  and 
parallel  operations  to  the  network,  one  obtains 
1 


A  =  (1  +  A,) 


—  (25) 

Cs 


where 


(2m2  +  lnnr)s  +  (8m2  +  8mr) 

Ax  = — (26) 

(m2  +  2mr)  s2  +  (Ipnr  +  Inn2  +  2m  +  2r)  s  +  8m 

Therefore 

[(m2   +   2mr)s2   +    (6m2   +   8mr  +   2m  +   2r) s   +    (8m2   +   8m  +   8mr)J 


A   = 


j_(m2C   +   2mrC)s3   +    (6m2C   +   8mrC   +   2mC   +   2rC   t  m2   +   2mr)  s2 
+    (8m2C   +   8mC   +   8mrC   +  lp-.r  +  Iuti2   +   2m  +   2r)  s   +   8m  J 

=—  (27) 

Da 

Similarly, 

B   =   1   +   B-l  (28) 

where 


19 


2ms  +  l_pn 
ms2  +  (2m  +  1) s 
Therefore 


BX  =  — ^ (29) 


ms2  +  (Ipn  +  1)3  +  ipn   Nb 

B  = =  —  (30) 

ms2  +  (2m  +  1) s      Db 

From  Eq.  (27)  and  Eq.  (30),  one  obtains 

B  -  A    Nb  .  Da  -  Na  .  Db 

=  (3D 

B  +  A    Nb  •  Da  +  Na  •  Db 
where 

Nb  ■  Da  =  (m3C  +  2m2rC)s^  +  (lOm^c  +  l6m2rC  +  3m2C 

+  i;mrC  +  re?   +  2m2r) s^  +  (36m3c  +  22m2C  +  2mC 
+  l±Qm   rC  +  l6mrC  +  2rC  +  &n?   +   3m2  +  12m2r 
+  L|.mr)s3  +  (56m3C  +  lj.8m2C  +  8mC  +  61^m2rC 
+  l6mrC  +  20m5  +  20m2  +  2l;m2r  +  12mr  +  2m 
+  2r)s2  +  (32m5C  +  32m2C  +  32m2rC  +  l6m3 
+  l^Om2  +  l6m2r  +  8mr  +  8m)  s  +  32m2      (32) 

Na  •  Db  =  (m3  +  2m2r)si+  +  ( 8m3  +  3m2  +  12m2r  +  1+jnr)  s3 
+  (20m3  +  18m2  +  2lpn2r  +  12mr  +  2m  +  2r)  s2 
+  (l6m3  +  2lpi2  +  l6m2r  +  8mr  +  8m)  s      (33) 

Bartlett's  representation  theorem  for  a  symmetric  network 
yields  the  transfer  function 

B  -  A    Nb  -  Da  -  Na  •  Db    N(s) 

T(s)  =  =  =  (3W 

B  +  A    Nb  .  Da  +  Na  .  Db    D(s) 
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N(a)  =  (m^C  +  2m2rC)s5  +  (I0m3c  +  l6m2rC  +  3m2C 

+  l|mrC)s^  +  (36m3c  +  22m2C  +  l4.8m2rC  +  l6mrC 
+  2raC  +  2rC)s-3  +  (56m5C  +  l|8m2C  +  6l4in2rC  +  8inC 
+  l6mrC  +  2m2) s2  +  (32m3C  +  32m2C  +  32m2rC 
+  l6m2)s  +  32m2  (35) 

D(s)  ■  (m3c  +  2m2rC)s^  +  (10m3c  +  l6m2rC  +  3m2C  +  lj.mr 
+  2m3  +  l|m2r)s^  +  (36m3c  +  22m2C  +  U8m2rC 
+  2l|m2r  +  l6m3  +  6m2  +  8mr  +  l6mrC  +  2mC 
+  2rC)s3  +  (56m3c  +  lj.8m2C  +  6luu2rC  +  8mC  +  l6mrC 
+  l|0m3  +  38m2  +  I|_8m2r  +  2l\mr  +  l\m  +   [|.r)s2 
+  (32m3c  +  32m2C  +  32m2rC  +  32m3  +  6i|m2  +  32m2r 
+  l6mr  +  16m) s  +  32m2  (36) 

In  order  for  this  transfer  function  T(s)  to  have  a  notch 
in  its  frequency  response,  the  values  of  m,  r,  and  C  must  be  so 
chosen  that  the  numerator  polynomial  N(s)  has  a  factor  of  the 
form  ( s2  +  to2).   The  values  of  m,  r,  and  C  can  be  determined  by 
dividing  N(s)  into  even  and  odd  polynomials  which  are  subjected 
to  a  Routh  array  calculation.   First  of  all,  let  N(s)  be  de- 
noted as: 

N(a)  =  a^s^  +  a^s^  +  8333  +  S2S2  +  8]_s  +  &q  (37) 

The  Routh  array  yields 
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s3 


a5 


^3ak   "  82a5 
ai| 

pla2  "  ftl«l* 


=  P-> 


a3 

a2 

aia^  -  a0a£ 


=  Qi 


=  Ri 


RlQl  "  pla0 

0  0 

There  must  be  a  row  of  zeros  at  the  end  of  this  Routh  array  in 
order  to  find  a  common  factor  in  these  polynomials,  namely, 


(38) 


Substitution  of  these  identical  representations  of  P-i,  ft-.,  and 

Rx  into  Eq.  (38)  yields 


(3384  -  a2a^)  |_(ai«4  -  a0a^)a2  -  (838^  -  a2a^)a0| 


=  (ana),  -  ana 


l8^  -  a0a5;   •  a^ 


(39) 


The  other  factor  which  determines  the  notch  frequency  00  is 


Ql       0 
(a2  -  —  •  Bu) s^  +  a0  =  0 


(ko) 


[a2(a3alj.   "    a2a^)    -    a^a-^   -    a0a5) 
+    ag^a^   -    a2a^)    =   0 


(111) 


On  comparing  coefficients  of  polynomials  in  Eqs.  (35)  and  (37), 
one  identifies 
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a^  =  rn^C   +   2m2rC 

a^  =   lOm^C   +   3m2C   +   l6m2rC   +  IjmrC 

83   =   36m3c   +   22m2C   +  l|.8m2rC   +   l6mrC   +   2mC   +   2rC 

a2   =   £6m3C   +  J+8m2C   +  61|in2rC   +   l6mrC   +   8mG   +   2m2 

a1   =  32m3c   +   32m2C   +   32m2rC   +   l6m2 

a0    =   32m2 

In  order  to  achieve  Eq.  (39),  several  steps  are  calculated 
separately  as  follows: 

(aialj.   "    80aS)     •    a2 

=    [(32m3c   +   32m2C   +   32m2rC  +   16m2)     ■    ( lOm^C   +   3m2C   +   l6m2rC 

+  IpirC)    -    32m2    •    (m3c   +  2m2rC)J    •     [(56ra3c   +  l|.8m2C   +   614m2 rC 

+   l6mrC   +   8mC   +  2m2) J 

=  32m3c    •    ["(560C2m6  +   1208c2m^  +   2096c2m5r  +   872C2m^ 
+   3528C2m^r  +   2560C2m1J-r2   +   21].8C2m3   +   1936C2m3r 
+   2880C2m3r2   +   1021j.C2m3r3   +   2i4.C2m2   +  ^2!|C2m2r  +  91|l4-C2m2r2 
+   5l2C2m2r3   +   32C2mr  +   61j.C2mr2   +   6l4.C2mr3)    +    ( 2^Cm5 
+   302CiA  +   61j4CliA-r  +   110Cm3   +   606Cm3r  +  i|l6Cm3r2   +   12Cm2 
+   176Cm2r  +   232Cm2r2   +   l6Cmr  +  32Cmr2   +   32Crar3) 
+    (8m^  +  3m3   +   12m3r  +  i(.m2r)  J  (lj.2) 

(a3ai+   "    a2a5>     '    a0 

=    [(36m3c   +   22m2C   +  l|8m2rC   +   l6mrC   +   2mC   +   2rC)     •    ( 10m3c 

+  3m2C   +   l6m2rC   +  IjmrC)    -    (£6m3c   +  lj.8m2C   +   62+m2rC   +   l6mrC 

+   8mC   +   2m2)     •    (m3c   +   2m2rC)]    •    32m2 

=   32m3c    (30l|Cm5  +   2800:114  +   SSOCm^r  +   688Cm3r  +   78Cm3 
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+   6l+OCm3r2   +   6Cm2   +   172Cm2r  +  l+l6Cm2r2   +   ll+Cmr  +   96Cmr2 
+   8Cr2    -    2iA   -   l+m3r)  (1+3) 

'aial+  "  aoa5'a2  "   'a3al|  "  a2a5'ao 

=  32m3c    [(560C2m6   +   1208c2m£  +   2096c2m5r   +   872C2m14-  +   3528c2m^r 
+   2^bOG2m^r2   +   2l).8c2m3   +   1936C2m3r  +   2880C2m3r2   +   102l+C2m3r3 
+   21+C2m2   +  l|.21+C2ra2r  +   9l+l+C2m2r2   +   5l2C2ra2r3   +   32C2mr 
+  61+C2mr2  +  6l+C2mr3)    +    ( -60Cm^  +   22CnA   -   236CreAr  +   32Cm3 

-  82Cm3r    -    22l+Cra3r2   +   6Cra2    +  1+Cm2r    -    l81+Cm2r2    +    2l+Cmr 

-  61+Cmr2   +  32Cmr3   -    8Cr2)    +    ( 10m^  +   3m2   +   l6m3r 

+  l+m2r)]  (1+4). 

(aia^   -    a0a^)       •    a^ 

=   61+m^C2    •    (20Cm^  +   26Cm3   +   52Cm3r  +   6Cm2   +  l+6Cm2r 
+   32Cm2r2   +   8Crar  +   8Cmr2   +   8m3   +  3m2   +   12m2r  +  1+mr)2 
•    (l+OCm3   +   61+Cm2r  +   12Cm2   +   l6Cmr)  (1+5) 

Multiplying  Eq.    (1+1+)    by    (aoar    -    &2a$    an<^    expanding   it,    then 
combining  with  Eq.    (1+5)  ,    one   obtains   a    polynomial   equation 
H(m,    r,    C)    identical    to  Eq.    (39)    such   that 

H(m,    r,    C)    =  ^C3   +   (3C2+    /  C   +   6    =  0  (1+6) 

where 

-<  =  69120m11  +   215616m10  +  l+56l92m10r  +   27l+383m9 

+   127l+176m9r  +   1198080m9r2   +   185096m8  +   ll+51+096ra8r 
+   2929920m8r2   +   l561+672m8r3  +   70816m?  +   8828881^ 
+   2932096m7r2   +   321+6o80m7i-3   +   105l808mV<-  +   15216m6 
+   306880m6r  +   l559888m6r2   +   2758ll+l+m6r3   +   1703936m6r^ 
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+  262lJ4l|.m6i-5  +    I680m5   +   60781pn5r   +   l4-77320m5r2   +    1233280m5r3 

+  Il67360m5r^  +   327680m^r5  +   72iA  +   6320m^r  +   8l|432m^r2 

+  313920^1^3   +  1;0960cWJ-i4  +   l638[|Wj-r5  +   26l4m3r  +   8056m3r2 

+  I(.6032m3r3  +  77l81|m3r^  +  l|0960m3r5  +  320m2r2   +  3680m2r3 

+  7296m2i4  +  5l20m2r5  +   128mr3   +   256mA  +   2$bmr5 .  (lj.7) 

p  =   -22Ij.80m10   -   3151|V9    -    1381|ij.8m9r   -    11780m8   -    2102l6m8r 

-  3l8208m8r2   +   1358m7    -    8[|.196m7r   -   l|2!j.928m7r2   -   323582pn7r3 
+   1698m6   -    8090m6r   -    20791Utm6r2    -   !|078o8m6r3   -   12880m6r1+ 

+  330m^  +  3352m^r   -   ij4088m5r2   -   206208nA-3    .    133376^^ 
+   l8m^  +   830r>Ar   -   3206m'+r2   -   59808m^r3   .   69888m1+r1+ 
+   1021^4-^  +  l).8m3r  +   112m3r2   -   I5568m3r3   -    23872m3r1J- 
+   6656m3r5  +   l^m2r2   -    102l+m2r3   -    52ij.8m2r^  +    l536m2r^ 

-  I|8mr3    -   6i|0mi4  +   128mr5   -   3  2  A  (i|8) 

/=   -980m9    -   85tpi8   -   l).588m8r   -   158m7    -   i^Wr   -    7H2m7r2 
+  33mb   -   928m6r   -   6692m6r2   -   36l4-8m6r-3   +   c,m5   .   g^5P 

-  2552m^r2   -   350lim5r3   +   33^  +  68m^r2    -   832m^r3    _   6^1^ 
+  l|.0m3r2   +   l6m3r3   +   l6m2r  (lj.9) 

6   =   -10m8   -   3m7   -   36m7r   -   10m6r   -   32m6r2   -    8m^r2  (50) 

Returning   to  Eq.     (J4.I)  ,    one    identifies 

82(838^   -    a2a£) 

=  [(56m3c   +  lj.8m2C   +   64ra2rC   +   l6mrC   +   8mC   +   2m2)] 

•    [2mC(l52m5c   +   ll^Om^C   +  I^OiArC   +   ii^m^rC   +  39m3c 

+   320m3r2C   +   3m2C   +   86m2rC   +   208m2r2C   +   7mrC   +  lj.8mr2C 

+  l|r2C    -   wh   -   2m3r)J  (51) 
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=   2mC-(l60m5c   +   208nAc   +  l|l6m^rC   +  I|.8m3c   +   368m3rC 
+   256m^r2C   +   61j.m2rC   +   6[j.m2r2C   +   6L±nM  +   96m3r   +   2hp? 
+  32m2r)    •     (I0m3c   +   l6m2rC   +   3m2C   +   IpnrC)  (52) 

ao(a3al4.   "    a2a5> 

=   l|m2C    •     (2lj.32m6C    +    22l4.0m5c    +    70lj.Om5rC    +   550lpArC    +   621^^ 
+   5l20m^r2C   +  lj.8m3c   +   1376m3rC   +   3328m3r2C   +   112m2rC 
+   768m2r2C   +   6l4mr2C   -    l6m^   -   32nAr)  (53) 

Substituting  Eqs.     (5D  ,    (52),    and    (53)    into  Eq.     ( J4.I)    and 

o 

expanding   them   term  by   term,    one    obtains   W    . 

P]_C   +   qx 

co2   =  r± (5W 

P2C^    +    tgC    +    q2 

where 

px   =   21|32m6   +   22l|0m^  +   70lj.0m^r  +   550lpiAr  +   621|jiA  +   5l20ml*-r2 

+  l|8m3   +   1376m3r  +   3328m3r2   +   112m2r  +   768ra2r2   +   61|mr2 

(55) 

qx  =   -l6m^   -   32m^r  (56) 

P2   =  3^56m7   +   6288m6   +   1382i4Bi6r  +  I|508m5  +   21l\l\.Om^r 

+   l8[|.32m5r2   +   l508rtA  +   13120irAr  +   23552raJ+r2   +   8l92m^r3 
+   228m3  +  3852m3r  +   11776m3r2  +   8l92m3r3   +   12m2   +  536m2r 
+   2880m2r2   +   3072m2r3   +   28mr   +   31+i+mr2   +   5l2rnr3   +   l6r2 
+  32r3  (57) 

t2   =   -196m6   -   lOOra^   -   6l4.8m^r   -  vh  -   336m^r   -    5l2ml+r2   +   3m3 
-    I8m3r   -   256m3r2   +   7m2r   -    l6m2r2   +  Ijinr2  (58) 


26 


q2 


-rri' 


5   -    2A-  (59) 


By  the  argument  used  previously,  one  raust  have  the  maximum 
value  of  C  in  Eq.  (I4.6)  .  It  is  clear  that  one  can  take  the  par- 
tial  derivative   of  H(m,    r,    C)    with  respect   to  m   and   obtain 

3E  20  ,   2+  dG  d$  3Z 

—  =   3^C^   —   +   G^  —   +    2(3C   —  +   CT   —  +    7  — 
dm  3m  %m  dm  gm  dm 

3/       $$ 
+  C  —  +  —  =  0  (60) 

»m       3m 

dG 
Putting  —  =  0   in  Eq.    (60)    yields 
0m 

d*      ,        3P      ?        3?  25 

—  CJ   +   —   Gd    +  —  C+  —  =0  (61) 
dm              dm              3m            3m 

In  order  to  simplify  the  notations,  some  symbols  used  in 
the  following  pages  are: 


xl 

8+ 

3m 

jll 

•<11 

3  m2 

S   J. 

■*12 

2mdr 

h 

2?> 

1m 

Pll 

2?2P 
3  m2 

Pl2 

22P 
3m0v 

"2 

9r 

-€± 

^22 

??2 

SG 

Cl 

2m 

P2 

3? 
3v 

P22 

2»2(3 
3r2 

3G 

c2 

3r 
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A 

'11 

32y 

<?m2 

^12 

.Pm^r 

,96 

5] 

3  m 

*25 

611 

^m2 

512 

^26 

oy 

y  = 

x2 

.9? 

22 

_£V 

<?r2 

^2C 

C12 

<?m?r 

?6 

52 

?r 

•?26 

522 

=  ^2 

<?2C 

Cll 

#ra2 

Returning  to  Eq.  (1+.6),  H(m,  r,  C)  = -<c3  +  (3C2  +  /C  +  6  =  0, 
one  can  see  that  this  is  a  cubic  equation  of  C  in  terms  of  m 
and  r.   Generally,  a  cubic  equation  can  be  solved  using  Car- 
dan's formulas  (11)  which  are  described  in  the  literature. 
However,  all  coefficients  of  C  in  Eq.  (I4.6)  ,  obviously,  are  com- 
plicated polynomials  in  terms  of  m,  r.   Thus  it  would  be  diffi- 
cult to  solve  for  C  in  terms  of  m,  r  by  the  classical  method. 
Also,  since  it  is  desirable  to  use  the  computer  for  calculation, 
the  complexity  of  the  expression  for  C  would  lead  to  problems 
with  computer  storage  capacity. 

To  avoid  these  problems,  an  optimum  process  based  on  the 
gradient  vector  method  (Ref.  9)  is  applied  to  solve  for  the 
specific  values  of  m,  r,  and  C  in  Eq.  (l|i>),  in  which  C  is  a 
maximum  value.   This  method  is  described  below. 

Rearranging  Eq.  (ij.6)  yields 


m. 
arms 
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-    -(/C  +  5) 
C2  =  (62) 

-<C  +  p 

In  order  to  obtain  the  maximum  value  of  C,  one  must  take  the 

first  and  second  partial  derivative  of  C  with  respect  to 

Then  putting  C^  =  C-q  =  0,  yields  a  linear  form  of  C  in  tei 

of  m,  r.   Several  steps  are  achieved  as  follows.   First  of  all, 

taking  the  partial  derivative  of  C  in  Eq.  (62)  with  respect  to 

m  yields 

UC   +  PM-ZtC  -  "/Ci  -  Si)  +  ("/c  +  b)(M1   +  ^]_C  +  P]_) 

20^  = 

UC  +  p)2 

(63) 

For  convenience  in  calculation,  one  puts  C,  =  0  in  Eq. 
(63)  before  proceeding  to  calculate  C-,-,.   Therefore,  Eq.  (63) 
becomes 

UC  +  P)(-^C  -  61)  +  (/C  +  5)UXC   +  (3t_)  ■  0  (6l|) 

Rearranging  Eq.    ( 6/4.)    yields 

-<]6C    +    p-i/C    -    p/iC    -    -<6t_C    -    P6t_    +    p]_6 
r«£    —    


•*/l      "    -<-]*'' 


(65) 


Secondly,    taking   the   partial   derivative   of  C    in  Eq.    (65)    with 
respect   to  m  yields 

!.(■</]_    -   ^1/)(-<116C    +   -^SC-l   +    Pn/C    +    Pi/Ct_    -    P-/11C 

-  p/iCi    -   -<6l:lC    -    p5lx   +    Pn6)     -    Ut_6C    +    Pj/C 

-  p^C    -    ^6x0    -    P6t_   +    PxSjt^/n    -    -<1]L/) 

2CCX   =  , 

1  (66) 

Simplifying  Eq.    (66)    and   putting   C-j_   =  0   again,    one    obtains 
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[-<P?X151    -    -<P/x5ll   +   ^Pn^iB    -    ■<11^51   +   ^1(iY61 

-  -<Piyil6    +   -^llPl^o    -    ^lPll-Zs] 

C   = — 

[-<2/x6xx    -   -<2/xi6l   +   ^n^Bi    -    ^xx^xS    +   -^lP^ 

-  ^u/fi  +  -'iPii7'2  -  ^nPi^2  +  ^Pi^ii  -  ^lP^ii] 

A    +   B 

=    (67) 

D   +   E 

Similarly,  taking  the  derivative  of  Eq.  (62)  with  respect  to  r 

and  using  the  procedures  described  above  yields 

I  -<P7,22S2  "  -<P>/2522  +  "<P22?'25  "  -<22Pr'52  +  •<2P,/52 

-  -<P2')/225    +   -<22P2'>/5    -    -^2P22^J 

C  ■  ■» — 

|_-<2/2522    "   -<2/2252   +   -<-<22'/52    -   •<-<22'>/25    +   -<22P7^2 

-  -<P22r'y2    +   •<2P227'2    "    -<22P2'?/2   +   -<P2/'r22    -    -<2p7722J 

U   +   V 

■=     (68) 

P  +  Q 

Equating  Eq.  (67)  and  Eq.  (68),  one  obtains  a  polynomial  equa- 
tion G(m,  r)  in  terms  of  m  and  r,  i.e., 

A  +  B    U  +  V 

G(m,  r)  = =  0  (69) 

D  +  E   P  +  Q 


A  +  B 

Substitution  of  C  =  into  Eq.  (6l)  yields 

D  +  E 


A  +  B        A  +  B         A  +  B 

Htm,  r)  =  -<x( )i   +  Px( )   +  7i( )  +  61  =  0   (70) 

D  +  E        D+E         D  +  E 
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The  value  of  C-jn  (i.e., )  shall  be  investigated  to 

d  mar 
assure  that  C  has  a  relative  maximum.   If  C-.p  =  0>  then  C  has 

a  saddle  point.   If  C-jp  ¥   °  (greater  or  less  than  zero),  then 

C  has  a  relative  maximum  or  minimum. 

Returning  to  Eq.  (63), 

UC  +  p)(-7<C  -  Ycx    -  St)  +  (/C  +  eJUCi  +  rf-iC  +  Px) 

2CC   =  = 

1  UC  +  p)2 

(71) 

Taking   the    partial   derivative    of  Eq.     (71)    with  respect    to   r 
yields 


2(CC,9   -   C-,C9) 


—  rue  +  p)(-7ic  -  rtj  -  Oifl 

0r    L  -1 


'12        °lu2;  ,    „         „,2 

(-<c  +  pr 


UC   +   p)(-7]_C   -  /cx  -   6X)    •   —   Uc  +  p)2 

9r 

UC    +    p)^ 

3 


[Vc    +   6)(^C!    +   rfjC   +    px)] 


(-tO  +  p)2 

(/C    +    6X^0!    +   ^]_C    +    px)     •    —    UC    +    p) 

3v 


UC  +  p)^ 

(72) 


Expanding   and   rearranging  Eq.    (72),    one    obtains 


•<6    -    p/ 


2CC-.0    "  o    '    C-|p 

(-<C  +  p)2         12 
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C    (-<i2/?/  "   ^\2    +  J-\^2   +   -<2'ri)    +   C(-<152   +   "<251 

UC    +    p)2    L 

+   ^2^1  +   Pl^2   +   •<125    "   "<512   +   Pl2^"    P*12'    +    (Pl52  +   P25l 

-,                   2 
+    (3126    -    (3612) r    C/O   +   SjUjC    +    P^  U2C   +    ?2>  ' 

(73) 


Thus 


c12  =  — -g J  |c^(^12y  -  -<y12  +  ^Y2  +  ^2ri) 

2c(^c  +  p)    -  ^s  +  pr  [L 

+   C(^152    +   -<261   +    P2^1+    Pl?2   +   "<125    "    "<512    +    Pl2^ 

-  pyi2)    +    (px52    +    p251    +    pi26    -    P612)J 

-    (Yc  +  s)  Uuc  +  Pi)  (^2C  +  Pa)  I  I 

L-<C    +   p  JJ 

WX   -   WY 
= (7W 

wz 


Among  those  solutions  of  C,  m,  r,  C12,  in  Program  1,  one 
educated  guess  is 

C  =  9185.78979  x  10"6 

—  9185.79  x  10~6   farad 
m  =  0.30l|2 
r  =  0.0179   ohm 
C12  =  ^7«i+31   farad  per  ohm. 
For  this  set  of  solutions,  C-,2  Is  greater  than  zero,  so  the 
dependent  C  possesses  a  maximum. 

Substitution  of  these  values  of  C,  m,  and  r  into  Eq.  (3i+) 
yields 
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N(s) 
T(s)    = 


D(s) 

where 

N(s)   =   (0.030537^2s^  +  0.59999283^  +  3.837814.82s3 

+   28.7633238s2   +   I65.081;l871s   +  322.3679696)    x  C 

D(s)   =  (0. 03053714.2  s5  +  7.14.5031213^  +  122.38l4.06lj.is3 

+   677.1501976s2   +   1291.8172755s   +   322.3679696)    x  C 

Furthermore,    substituting   s   =    jco  yields 

[(0.5999928oA  -   28.7633238co2  +  322.3679696) 
+  j (o.03o537l|.2a)5  -  3.837814.82a)3  +  165.08I4.18710)] 


T(jto)    = 


[(7.14.5031210)^  -  677.15019760)2  +  322.3679696) 

+  j (0.03053714-20)5  -   122.38I4.06I4.I0)3  +  1291.8l72755u)_l 


[(0. 5999928a)11  -   28.76332380)2  +  322.3679696)  2 

.2         +   (0.03053714-20)5  -  3.837814.82a)3  +  l65.o8l4.i87lo))2J 
T(jco)        =  ■= : 5 5 

[( 7. 14.5031210)'+  -  677.1501976a/  +  322.3679696)^ 
+  (0.03053714-20)5  -  122.38l4.06l4.ia)3  +  i29l.8i72755o))2J 

By  inspection,    one   can    see      that 

i  2 

T(  jo))  j      =   1  at      o)  =  0 

|t(  jo>)  I      =   1  at      to  =  °^ 

The  log-log  plot  of  T(jto)    versus  to  is  shown  in  Fig.  12. 
This  frequency  response  shows  that  the  transfer  function  T(s) 
represents  a  notch  network  with  equal  amplitudes  at  zero  and 
infinite  frequencies.   Data  for  this  plot  are  obtained  with 
Program  I4.  in  Appendix  C. 
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Pig.    9.      An   equal   amplitude   notch  filter 
resistance-capacitance   network. 
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Fig.  10.   Mid  short-circuited  impedance  in  Pig. 
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Fig.  11.   Mid  open-circuited  impedance  in  Fig.  9. 
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COMPARISON  OP  NOTCH  FILTER  Q'S 

The  parameter  Q  represents  the  sharpness  of  the  circuit 
frequency  response  curve.   It  can  be  defined  as 

energy  stored  in  circuit 

Q  =  2jt  .  — . 

energy  dissipated  in  circuit  during  one  cycle 

Similarly,  the  factor  Q  defines  the  sharpness  of  the  notch 
spectrum.   The  bandwidth  of  a  notch  filter  is  commonly  defined 
as  the  width  of  the  band  of  frequencies  over  which  the  output 
power  does  not  drop  to  less  than  one-half,  or  -3  db  of  the  out- 
put power  at  notch  frequency.   The  explicit  formula  is 


r0  u0 

Q  = 


<Afo>3db   (A"o>3db 


where  fQ  is  the  notch  frequency  and  Af0  is  the  number  of  cycles 
off  notch  frequency  at  which  the  responses  are  70.7  per  cent  of 
its  peak  value.   In  other  words,  this  factor  Q  is  the  reciprocal 
of  the  bandwidth  in  units  of  the  notch  frequency. 

Case  1.   Equal  amplitudes  notch  filter  of  the  Goldman  type. 

From  Fig.  i|,  which  is  the  response  for  equal  amplitudes 
notch  filter  of  the  Goldman  type,  one  obtains 

1.55 

Qn   =  =  0.0803 

19.55-0.192 
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Case  2.   Unequal  amplitudes  notch  filter. 

The  bandwidth  is  undefined  because  of  the  unequal  magni- 
tudes of  the  curve  (Pig.  8)  at  zero  and  infinite  frequencies. 
Therefore  Q2  for  this  filter  is  undeterminable. 

Case  3.   Equal  amplitudes  notch  filter. 

6.82 

0.0193 


3   352.3  -  0.18 


Qo  is  smaller  than  Q, ;  the  notch  for  the  simple  Goldman 
type  network  is  sharper  than  that  of  the  complicated  one.   The 
smaller  Q-i  for  the  complicated  network  indicates  that  a  wider 
bandwidth  is  achieved;  the  bandwidth  is  increased  by  Jb   per 
cent.   Increased  bandwidth  is  desirable  in  carrier-frequency 
servomechanisms  subjected  to  wide  bandwidth  input  signals. 
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CONCLUSION 

An  RC  notch  filter  has  a  zero  minimum  frequency  response 
curve.   Three  cases  have  been  exhibited  in  this  paper: 

1.  An  equal  amplitude  RC  network  of  the  Goldman  type. 

2.  An  unequal  amplitude  RC  network. 

3.  An  equal  amplitude  RC  network  (see  Fig.  9). 
Determination  of  proper  parameters  for  the  existence  of  a 

notch  frequency  in  each  case  is  emphasized.   Particularly  in 
the  third  case,  an  optimum  process  is  applied  to  decide  the 
specific  values  of  m,  r,  and  C  for  the  notch  frequency  to  occur 
at  6.82  radians  per  second.   (See  Program  1.) 

The  graph  for  the  frequency  response  in  each  case  shows 
its  particular  amplitude  characteristic. 

Generalized  Goldman  RC  notch  networks  have  led  to  a  maximum 
problem  of  an  unorthodox  type.   Specific  procedures  have  been 
exhibited  for  reducing  this  maximum  problem  to  the  problem  of 
solving  two  bivariate  polynomials.   Numerical  trial  and  error 
procedures  affecting  numerical  solutions  on  a  digital  computer 
are  given  in  Appendix  C. 
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APPENDICES 
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APPENDIX  A 


Description  of  the  IBM-1620  Computer  Program  Used  for 
Solving  Two  Bivariate  Polynomial  Equations 
in  Program  1. 


Several  numerical  methods  are  available  for  solving  non- 
linear polynomial  equations.   The  procedure  selected  was  the 
numerical  trial  and  error  method. 

From  Eqs.  (69),  (70)  in  this  paper,  one  obtains 

A  +  B    U  +  V 
G(m,r)  = =  0  (75) 

D  +  E    P  +  Q, 

A  +  B  ,       A  +  B 

H(m,r)  =  -<-,( )J  +  Sn( r 

D  +  E        D  +  E 

A+B 
+  7M )  +  51   =  0  (76) 

D  +  E 

where  all  parameters  in  Eqs.  (75)  and  (76)  are  high  power  poly- 
nomials in  terms  of  m,  r. 

First  of  all,  let  some  symbols  used  in  Program  1  be  de- 
noted as: 

AX  =  ^  CX  =  Y 


AMI  =  ax   =  —  CM1  =  \   =  — 

Jm  ^m 

a2*  d2r 

AM2  =  .<, ,  =  CM2  =  /, ,  =  

11   ?m2  1X   dm2 

Z2^  ?** 

AM1R1  =  -<12  =  CM1R1  =  YX2   =  - 


$ra0r  2m3v 


AR1 

=   j.2    = 

3-c 
9* 

AR2 

=  -c22 

32* 
"  gT2 

BX   = 

=   P 

BM1 

=   Pi  = 

2(3 

#m 

BM2 

"   011 

2 
2^B 

?m2 

52P 

BM1R1   =    p12   =  

3m3r 

BR1 

=   P2   - 

3P 

3r 

BR2 

"    P22 

22P 

DM  = 

!  A  m 

CR1 

=  Y2  =  — 

3v 

CR2 

_  Y       —  

-   '22   ~        p 

DX 


DM1   = 

51   =  ~ 

dm 

d26 

DM2    = 

5ll   =  — 

<?m 

a25 

DM1R1 

—  ^     — 

-  o12  - 

35 

DR1   = 

52  =  t; 

2r 

3^5 


DR2 


;22  - 


3^ 


DR   = 

Ar 

C12    : 

=  c12 

= 

32C 

Smdv 

GX    = 

G(m, 

r 

) 

HX   = 

Htm, 

r 

) 

k3 


XM  =  m 
R  =  r 

Because  Eqs.  (75)  and  (76)  contain  high  power  polynomial 
functions  of  m  and  r,  an  extremely  large  computer  capacity  is 
required.  Since  the  maximum  storage  capacity  of  the  IBM-1620 
using  Fortran  II  is  limited  to  60000,  the  program  for  solving 
the  equations  was  separated  into  two  subprograms  (shown  in 
Program  1) .  ■ 


kk 


The  numerical  trial  and  error  method  is  illustrated  as 
follows. 

Let  (rag,  Tq)    be  an  initial  approximation  to  a  root  of  Eqs. 

(75)  and  (76),  and  let  Am  and  Ar  be  the  increments  of  m,  r 
respectively.   Substitutions  of  m0,  r0  into  Eqs.  (75)  and 

(76)  yields 

G(m0,  r0)  =  g(m0,  r0)  (77) 

H(m0,  r0)  =  h(m0,  r0)  (78) 

Theoretically,  both  function  g(m0,  r0)  and  function  h(m0,  r0) 
should  be  equal  to  zero  if  (mn,  r0)  is  an  exact  root  of 
G(ra,  r)  and  H(m,  r) .   This  is  difficult  to  achieve  because  the 
convergent  point  is  rather  hard  to  find  for  high  degree  poly- 
nomial equations.   Therefore  the  increments  Am  and  Ar  are 
added  to  make  G(m,  r)  and  Htm,  r)  possibly  convergent.   Further- 
more, other  initial  approximations  should  be  tried  to  get  all 
possible  convergent  points.   Care  must  be  taken  with  the  roots 
of  G(m,  r)  and  H(m,  r) ,  since  there  can  be  both  real  and  complex 
ones.   However,  only  real  roots  are  of  interest  in  this  work. 

Referring  to  the  flow  chart  (Appendix  C)  for  Program  1, 
one  should  supply  both  initial  approximation  and  increments  of 
m,  r  and  set  the  iteration  =  N(l,  2,  3,  .  .  . ,  n) .   Then  one 
may  compute  polynomials  and  obtain  computer  results  in  sub- 
program A.   These  results  are  then  introduced  into  subprogram  B. 
Computation  of  these  polynomials  in  subprogram  B  yields  the 
required  solutions:   XM,  R,  GX,  HX,  C,  C12. 

There  will  be  11  roots  because  G(m,  r)  and  H(m,  r)  are  ll**1 
degree  polynomial  equations.   The  optimum  process  is  applied  to 
select  the  required  roots. 
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APPENDIX  B 

Description  of  the  IBM-1620  Computer  Program  Used 
for  Program  2,  Program  3,  and  Program  ly. 

IBM-1620  FORGO  is  applied  to  calculate  the  evaluation  of 

I      I  2 

|T(ju)|   versus  co  (angular  frequency)  for  the  equation  such 

that 

,      2    anun  +  an-l"n_1  +  •  •  •  +  alU  +  a0 

T(jco)!   = (79) 

bncon  +  bn.!^""1  +  .  .  .  +  b1co  +  b0 

Referring  to  Program  2,  Program  3,  and  Program  I4.,  one  can  easily 
obtain  the  evaluation  of  T(jco)    for  change  of  co  from  zero  to 
infinity.   Computer  results  are  shown  respectively  in  those 
programs. 


APPENDIX  C 


IBM-1620  Computer  Programs  Used  in  This  Work 


J+6 


(A) 


r 


READ  N 


[ 


0 


READ  INITIAL 
ESTIMATES  XM,  R 


f^READ  INCREMENT 
DM     PR  I 


< 


do  55 

1  =  1,' 


> 


XM  =  XM  +  DM 
R  =   R  +  DR 


COMPUTE 

POLYNOMIALS 


PUNCH:  XM  R 
AX  BX  CX  DX  AMI 
AM2  AM1R1  AR1 
.R2  BM1 


RETURN  TO  START 


(B) 


y 

READ: 

XM 

R 

AX 

BX  CX 

DX 

AMI 

AM  2 

AM1R1 

AR] 

AR2 

BM1 

COMPUTE 
POLYNOMIALS 


Flow  chart  for  calculation  of  two  bivariate 
polynomial  equations  in  Program  1 
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PROGRAM  1.  NUMERICAL  TRIAL  ■' '  FCR  CALCULATION 

OF  POLYNOMIAL  EQUATIONS 

C  C       CALCULATION  OF  TWO  PIvAPIATF  POLYNOMIAL  FGUATIONS 
*  1  2  0  5 

1  FORMAT (3E15. 5) 

21  FCRMAT(2( E15.5) ) 

2l    FCRMAT(5<  £15.5)  ) 

52  FORMAT! 13) 

5]  RFAD52»N 

5i     FCRMATI2F9.8) 

I  READ    ! »XM. 
RFADC(  l,C   '  ■ 
XM=XM-DM 
R=R-HR 

DO  551=1 »N 
X  M=XMH 
R=R+DR 

II  AX=69120.*XM**1 1+21 5616. *XM**10+456192.*XM**10*R+274384.*XM*«9+127 
14176. *XM**9*R+119c  .  -  «  *2  +  lb5096  .*XI'i**8  +  1454096.*XM**8*R  +  2 
29299  2i,.*XM**8*R**2  +  l  564672.  *XM**o*K«*  3  +  7081  6.  *XM*»  7  +  88  2888.  *XM»*  7* 
3R  +  2  932  v/96.*XM**7*R**2  + 3246080  .*XM**7*R»*  3  +  1  051 808  ;*XM»*7*R**4 

AX  =  Ax  +  15216.*XM**6  +  3  068  80.*XV.**6*K  +  15598  88.*XM**6*R**2  +  2  758144. 
l**6*R**3+17i  39  36.*XM**6*R**4+262144.*Xi-'**6«k**5  +  1680.*XM**5+60784. 
2*XM**5*R+477320.*XM**5*R**2+]  2  332  80.  *XM**5*R**3+  1 16  7360  .  *xK1**5«R** 
34+'-<2768  .  KXM*-*5*R**5+77.*XM-,t*4  +  63  2  0.*X''s*4*:<  +  R4437.SXM**4*R**? 

..      I  . "XM**4*R**3+4n9600.*XK**4»R**4+163640.*XM»*4*R**5 . 
l.*XM**3*R+8056.*XM**3*R**2+46032.  ;'  XM*  *3*R**3+77]  84.*XM**3*R**4H 
26C .*XM**3*R«*5+3  20.*XM**2*R**2+36  8o.*XM**2*R**3+7296.*XM»»2»R* 
3120.*XK**2«R**5+12  8.*XM*R**3+2  5  6.*XM*R**4+2  56.*XM*R**5 

BX  =  -22480.»XM**1 0-31  544. *XM**9-1 38  448.  *XM**9*R-tll  780. *XM»*8-2  1021  6 
1 .*XM**8*R-318208.*XM»*8*R**2+1358.*XM**7-84196.*XM**7*R-424928.*XM 
2**7^R**2-323584.*XM**7*R**3  +  1698.*XM**6-8090.*XM*tt6*R-207944.»XM»* 
36*R*»2-4U78o8.*XM**6*R**3-128  8  0.*XI-i**6*R**4  +  330.*XM**5 

hX  =  QX  +  33  52.*Xiv,**5*R-44088.*XM**5*R**2-2t,62  08.*XM**5*R**3-1333  76.*X 
1M»»5*R**4+18.*XM**4+830.*XM**4*R~3206.*XM**4*R**2-59808.»XM«*4*R** 
23-6  98  88.*XM**4*R**4+l 0  24 .*XM**4*R**5+48 . «XM**3*R+ ] 12.*XM**3*R*#2-] 
3  5  56T .*XM»*3*R**3-23872.*XM**3*R**4+6656.*XM**3*R**5+14.*XM«*2*R**2 

B  X  =  B  X  - 1 '  24.*XM**2#R**3-5248.*XM**2*R*«4+1  536.*XM**2*R*«5-48  .*XM*R* 
1*3-640.*XM*R**4+12  8.*XM*R**5-32.*R**4 

CX=-980.*XM*#9-854.*XM**8-4588.*XM**8*R-158.*XM**7-4194.»XM**7»R-7 
1112 .*XM**7*R**2+3  3.*XM**6-9  2  8 • *XM**6*R-66  92 .*XM**6*R*«2-3648 . *XM*» 
2  6*R**3+9.*XM**5-54.»XM**5*R-2  5  52.*XM**5*R**2-3  5  04.*XK«*5*R**3+3  3.* 
3XM-'  »4*R+6  8.*Xi'-,**4'*R)i*2-8  32.*XM**4*k*»3-64.*XM**4*R**4 

CX=CX+4o.*XM*<3*R**2+16.*XM»»3*R**3+16.*XM**2*R 

DX=-1I  .*XM**8-3.*XM**7-36.*XM**7*R-10.*XM#*6»R-32.#XM**6»R*»2-8.*X 
1K»*5*R**2 

AM]  =76'  321  .*X"**1  0  +  2  1  5616'  .*XM**9+4561920.*XM*#9*R+2469456.*XM#*8+ 
11 146  7584.*XM*#8*R+10  78  2  72  0.*XM**8*R**2+1480  7  68.*XM»»7+] 1632  76P. 
2*#7*R+2 34 39360. #XM»*7*R**2+1 251 7376. *XM**7*R**3+495712.»XM*»6+6180 
3216.*XM**6«R+7  0  5  246  7  7.*XM**6*R**2+22  722  56  0.*XM**6*R**3 
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AMI =AM 1+7 362  6^6. *XM**6*R**4+9 1296. *XM«*5+ 1841 280. *XM**5*R+9  3 59328. 

»*5*R**2+ 16548 8 64. *XM**5*R**3+ 1022361 6. *X  l**5*R**4+1572864.*XM* 

2*5*1  **5  +  84i.  .  ■      ■  ,   .   ,    ,  hUjoooui   a,.   h  .,   ^omosuu.-* 

'        •  4+1638400.  3+288. *XM**3 

L=AMl+25280.*XM**3*R+337728.*>       «*2+125ai    .  **3+1638 

1400.*XM**3*R**4+655360.*XM**3*R*       .    **2*R+24168. *XM**2*R**2+ 

2138o96.*XM**2*R**3  +  2  31<-   .  **4+122 88  0.*XM**2*R**5+640.*XM*R 

3**2+T36C.»XM*R**3+14592.*XM*R**4+]  240. *>  *R**5+128 .*R*«3 

AM]  =AM1  +  256.*R**4  +  ?56.*r' 

&M2=760320C.*XM**9+194      .  .  ■  **8*R+i 9755648 .' 

1*7+9174:  672. *XM**7*R+862617   .  =*2+lC36: 

2.*XM**6*R+1640  7552  0.*XM**6*R**2+87621632.*XM**6*R**3+2974272.*X 
+123148032.  -2+13633'     .    K-*5*R**3 

AM.'  =AM2+441 7  5936.  *XM**5*R**4+456480.*XM**4+92  06400.  *XM**4*R+467966 
140. *XM**4*R**;  +82744320.  <M**4         1 1 18080. *XK**4*R**4+7864320.* 
2XM**4*R**5+33600«*XM**3+1215680.  ■       4  ■  "346400. *Xi«'.**3*R**2  +  24665 
. *XM**3*R**3+23347200.*XM**3*k**4+63b36CC. *XM**3*R**5 

AM2=AM2+864.*XM**2+75B40.*XM**2*R+1013184.*XM**2*R**2+3767040,« 
1*2*R**3+49  15  2l  0.*XM**2*R»*4+196608O.*XM**2*R**5+1584.*XM*R+48336.* 
2XM*R**2+276192.*XM*R**3+463]   .  '      4+245760. *XM*R**5+640.»R**2+ 
37360. «R**3+1 4592. *R**4+1 0240. *R**5 

AMIR 1=4561 920. *XM**9+1 1 467 584. *XM**8+21 56 5  440. *XM**8#R+1 1632 768. *X 
*7+46B78 720. *XM**7*R+375 521 28. *XM**7*R**2+6 18021 6. *XM**6+4 104934 
24.*XM**6*R  +  6  81676oi.  .*XM**6*R**2+29450624.*XM**6*R**3+1841280.*X 
35+1871865  6.*XM**5*R+49646592.*XM**5*R**2+40894464.*XM**5*R*»3 

AM1R1=AM1R1+7864320**XM**5*R**4+3039   .  ■  - *4+4773200.*XM**4*R+184 
1992   .  ■  **4*R**2+23347200.*XM**4*k**3+8192    .    ;:  *4*R**4+25280  .  * 
2XM**3+675 456. *XM**3*R+3 767040 .*AM**3*R**2+6 55 3600. *XM**3*R**3+3276 
,*XM**3*R**4+792i*XM**2+48  336.*XM**2*R+414288.*XM**2*R**2 

AMlRl  =  AMlRl  +  9262ob.*Xi-**2*R**3  +  6144.  ,*XI  **2*  n»*4+1280 . *XK*K+22o80 
l.*XM*R**2  +  58  368.*X  M*R»*  3  +  51  2C0 .*XM*R**4+3  84.  *R*<-  2  +  1024.  *R«*3  +  1280. 
2*R**4 

AR1=456192.*XM**10+12741 76 . *XM**9+2 3  961 6C . *XM**9*R+ 1454096. *XM**8+ 
1585984'.,. *XM**8*R  +  4694016.*X^**S*R*1.  I  3  .«XM*#7+5fl64192.*XM»»7* 
2R+973824o.»XM**7»R*»2+4207232.*XM**7*R**?+306880.*XM»*6+3  119776.*X 
3M**6#R+8274432.»XM**6*R**2+6815  744. *X  -  -  '•  *  R**3+l  31  0  720.  *XM**6*R**4 

AR1.=  AR  1+60  784.  *XM**5 +  954640  .*XM**5*R  +  369  9  840. *XM**:L#R**2  +  466944  . 
iXM**5*R**3+163840o.*XM*1       ,  •     .  ■  **4+168864.*xM**4*R+941760. 
2*XM*.*4*R**2+16384o0.*XM**4*K**3+8  192ou'«*Xlvl**4*K**4+264.*XM**3  +  161  1 
32.*XM**3*R+138096.*XM**3*R**2+308   i  .    **3*R**3 

ARl=ARl+2C48u. .*XM**3*R**4+640.*XM**2*R+11040.*XM**2*R**2+29184.*X 
1M«»;  *R**3  +  2t>6>  0.*XM**2*R**4+384.*Xi"  *F  «  *2+1024.  *XM*R**3+1  280.*XM*R* 
2-4 

AR2  =  239616o.*XM**9  +  5  8  59  840.*XM«*a  +  9  388L32.*XM**8*R  +  5864  19  2.*X.M**7  + 
1 194764 8 i .*XM**7*R+ 12 621 696. *XM**7*R**2+31 19776. *XM**6+1 6548864. *XM 
2**6*R+20447;   .  ■  **6*R**2+524288u .*XM**6*R**3+954643.*XM»*5+73996 
,*XM**5*R+14     .  .  ■        »;         . 

AR2  =  AR2+168864.*XM**4+1  833[;'   .  .-  X  l**4*R#*2+32  76800 

l.*XM**4*R**3+16112.*XM**3+276192.*XM**3*R+92620l  .  •    3*R*»2  +  81' 
20.»XM**3*R**3+640.*XM**2+22080.*XM**2#R+87552.*>        **2+10240  . 
3XM**2*R**3  +  768.*X!':*R  +  30  72.*XM*R**2  +  312o.*XM*R**3 


U9 
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bMl=- 
11728. 
2496.* 
347664 
■■!--■ 
1*4*R* 
2*2-23 
3*2*R+ 

BM1 
1+3C72 
PUNCH 
PUNCH 
PUNCH 
CCNT I 
GO  TC 
END 


2248 
*XM* 

*XM 

Ml  +  1 
*3-6 

9232 
336. 

Ml  +  1 
I  x  ■' 
1»X 
2o 
2, 

NUb" 
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jC.*XM**9-2o3896.*XI"'*  '  .  >  XM**8*R-94240.*Xt-<**7-168 

-2  54  5664.*XM**7*R*  . *XM**6-589372.*XM**6»R-2974 

6»R**2-226r-u88.*/,  .     -'5-48540.  *XM**5*R-1 2 

**5*h!**2-244684S.  .  x  l**5*R**4 

,*XM**4+1«    .  <-22C44(  .»Xi»i"4»n»»2- 103104D«*Xn» 

.  ■  **4«3**44   .  ■     ,*Xf"'*3*K-12824.»XM**3*K-» 

*Xt^**3*H**3-279552.-,;>  •        R**5+144.*XM* 

*Xt-'**2*R**2-467C4.*X>'  ;  ;-71616.-'        :  *4 
9968.*XM**2*R**5  +  2o.*Xf/*1<s*2-2-48.-X-  i*R**3-10496.*XM*R**4 
*R*#5-48.*R**3-640.*R**4-128.*R**5 
M,R 
AXjBX.CX.DX.a  I 

.  ■   I1R]  »AR1  »AR2,   I 
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*1  20  5 

2C    FCRMAT(5( E15.5 ) ) 
1    FCRMAT(3E15.5) 
7    FCRMAT(2X,2HC=E15.5,5X,3HXI         '     .     , 5X .2HR=E1 5. 5 ,/ ) 

3  FCRMAT(2X,2HC  =  E15.5,5X,3HX'-  =  :-  '     .     -      1 ,4HC12=E12.  5  ,  /  ) 

4  FCRMAT(2X,2HR=F15.5  ,  6X  ,  ?HGx  =  c  1  ?  . '- •  C-X  ,  "J  Hr  IX  =  F  1  2  .  ;   ,/) 
K       "FAD       1,XM,R 

READ       2  UAX.BX.CX.PX.AM1 

RFAD      20iAM2»AMlRl»ARl  >AR2>I 

BM2=-2U23200.*XM**8-2271168.*XM**7-99682  56.  '■'•  X  ■'.**  7*R-6  5  9680.*XM**6- 
11 1772096. *XM**6*R-1 78 19648. *XM**6*R**2+b7036.*XM**5-3536232.*XM**5 
2*R-178469  76.*XM**5*R**2-13  59  0528.*XM**5*R**3+50940.*XM**4-242  700.* 
3XM**4*R-6  2  38  320.*Xf"'**4*R**2-122  34240.*XM**4*R**3 

dM2  =  BM2-36864lO.*Xivs*4*R**4  +  66l,0'*XM**3  +  67(;40.*XM**3*R-881760.*XM* 
1*3*R**2 -4 124 160. *XM**3*R**3-2 667520. *XM**3*R**4+216.*XM**2+99  60.*X 
2V**2*R-3  8  4  72.  -X •<**?* R** 2-7  1  7696  .  *xiv**2*R**3-8 3  86  56.  *XM*»2*R**4+12 2 
">PP.*X^**2*R**5  +  288.*XN'*R+6  7?.*Xf'*R**2-9  3  4  08.*XM*R**3 

I2=BM2 -1432 32. *XM*R**4+39936.*XM*R**5+28.*R**2-2048.*R»* 3-10496.* 
]  R**4  +  T-  72.*R**5 

1R1  =-12  46Ci3  2.*XM**8 -1681728.  *XM**7-5091 32  8.*XM**7*R- 589  372.*XM»* 
16-5  948  9  9  2.*XM**6*R-6  79  5264.*XM**6*K**2-4  8  54  0.*XM**5-2495  328.*XM**5 
2*R-7  34o54  4.*Xi';**5*R**2-294912  0.*X,-1**5*R**3-16  760.*XM**4-44u88  0.*X'>l 
3**  ,*R-3093120.*XM**4*R**2-2667520.*XM**4*R**3+3320.*XM**3 

BM1R1=BM1R 1-25648. *XM**3*R-71  76  96.  *XM**3*R**2-  11 182  08.  *XM**3*R**3+ 
24.*XM**2*R**3+99  840.*XM**2*R**4+5  6.*XM*R-6144.*XM*R**2-41984.*XM*R 
3**3+15  360.*XM*R**4-144.*R**2-2  5  6  0.*R**3+640.*R**4 

BRl=-138448.*XM**9-210?16.*X^«*8-6  364]6.*Xi-'**8»R-84196.*XM**7-r 
] 56.*XM**7*R-970752.*XM**7*R**2-8090.*XM**6-415888.*XM**6*R-1223424 
2.*XM**6*R**2-49152f'.*XM**6*R**3  +  3  3b2.*XI1'**5-88  176.*XM**5*R-618624. 
»*5*R**2-53  3  504.*XM**5*R**3+8  3u.*XM**4-6412.*XM**4*R 

BRl=BRl-1794  24.*XM**4*R**2-?79K52.*Xi-'**4*R**3  +  5120.*XM**4*R**4+4R. 
l*XiVi#*3  +  2  2  4.*XV!**3*K-46  704.*XN'**3*R**2-9b488.*XM**3*R**3  +  33280.*XM* 
2*3*R**4+28.*XM**2*R-3072.*XM**2*R**        ,J.*XM**2*R**3  +  7680.*XM**2 
~4-144.*XM*R**2-2  56  0.*XM*R**3  +  640.*Xf',*R**4-12  8.*R»*3 

uK2=-63  64  16.<  XM**8-8498 56. *XM**7-1941 504. *XM**7*R-4 15888. *XM**6-24 
146848. *XM**6»R-14745 60. *XM**6*R**2-B81 76. *XM**5-123 7248. *XM**5*R-1 
26C0512.*XM**5*R**2-6  412.*XVl**4-35  8b4a.*XM**4*R-8  386  56.*XM**4«R**2  + 
320481  .*X^**4*R**3+2  2  4.*XM**3-93  40  8.*XM**3*R-2864  64.*XM**3*R**2 

BR2  =  BR2  +  1331?i  .*Xf-'**  •3*R**3  +  ?8  .  *XM**2-6  1  44  .*XM**2*R-62976.*XM**2*R* 
1*2+30720. *XM**?*R**3-2  88.*XM*R-76  8  0.*XM*R**2+2  56  0.*XM*R**3-3  8  4.*R* 

1  =-8  8  2;,  ,*XM**8-6  83  2.*XM**7-36704.*XM**7*P-1  106  .  *XM*#6-293  58  .*XM* 

1  <6*R-497  84.*XN'**6*R**2  +  198.*XM**5-5568.*XM**5*R-40152.  *XM**5*R**2- 

221888. *XM**5*R»*  3  +  45  .*XM*«  4-2  70  .*XM**4*R-12  760.  *X'  I        2-17520.* 

3XM**4*R**3  +  13  2.*XN.*  *3*R  +  2  72.*Xi'-'**3*R**2-3328.*XM**3*R**3 

CM1=CM  1-2  56.  *XM**3*R**4+120.*Xl»l**2*R**  2  +  48.  *XM**2*R»*  3  +  32.  *XM*R**3 

CM2=-7u56  0.*XN**7-4  7824.*XM**6-2b6  92B.*XM**6*R-66  36.*AM**5-176146. 

l*XM**5*R-298  7u4.*XK**5*R** 2+990. *XM**4-2 7 840. *XM**4*R-2 0076 C.*XM»* 

24*R**2-10  9440.*XM**4*R**3+180.*XM**3-108C .*XM*#3*R-5 1040. *XM**3*R* 

3*2-7;  08.1  .  *  Xv;**3*R**3  +  396.*X^**2*K  +  ttl6.  »XM**2*K**2-9984.  *XM**2*R**3 

CM2=CM2-768.*XM**2*R**4+240.*XM*R**2+9  6.*XM*R**3+32.*R**3 


51 


CMlRl  =  -36  7«>4.*XM**7-29  358.*XM**6-99568.*XM**6*R-5  568.*XM**5-80304. 
1*XM**5*R-6  5664.*XM**5*R**2-27  0.*XM**4-2  5520.*XM**4*R-5  2560.*XM**4* 
2R**2+132.*XM**3+544.*XM**3*R-9984.*XM**3*R**2-1024.*XM**3*R**3+240 
3.*XM**2*R+144.*XM**2*R**2+96.*XM*R**2 

CR1=-4588.*XM**8-4194.*XM**7-14224.*XM**7*R-928.*XM**6-13384.*XM*» 
16*R-10944.*XM»*6*R**2-54.*XM**5-5104.*XM**5*R-10512.*XM**5*R**2+3  3 
2.*XM**4  +  136.*XM**4*R-2496.*XNt**4*R**2-256.*XM**4*R**3  +  8  0.*XM**3*R  + 
348.*XM**3*R**2+48.*XM**2*R**2 

CR2=-14224.*XM**7-13  384.*XM**6-21868.*XM**6*R-5104.*XM**5-21024.*X 
1M**5*R+136.*XM**4-4992.*XM**4*R-768.*XM**4*R**2+8  0.*XM**3+96.*XM»* 
23*R+96.*XM**2*R 

DM1=-80.*XM**7-21.*XM**6-2  52.*XM**6*R-60.*XM**5*R-192.*XM**5*R**2- 
140.*XM**4*R**2 

DM2=-560.*XM**6-126.*XM**5-1512.*XM**5*R-30  0.*XM**4*R-960.*XM**4*R 
1**2-1 60. *XM**3*R**2 

DM1R1=-252.*XM**6-60.*XM**5-384.*XM**5*R-80.*XM**4*R 

DR1=-36.*XM**7-10.*XM**6-64.*XM**6*R-16.*XM**5*R 

DR2=-64.*XM**6-16.*XM**5 

A=AX*BX*(CM2*DM1-CM1*DM2)+BM2*DX*( AX*CM1 -AM1*CX ) 

B=BX*CX*I  AM1*DM2-AM2*DM1)+BM1*DX*(  AM2*CX-AX*CM2 ) 

D=AX**2*(CM1*DM2-CM2*DM1)+AX*AM2*( CX*DM1-CM1*DX ) +AM2*BX*CX*CM1 

E=CX**2*< AM1*BM2-AM2*BM1)+CX*CM2*( AX*BM1-AM1*BX ) -AX*BM2*CX*CM1 

U=AX*BX*(CR2*DR1-CR1*DR2)+BR2*DX*( AX*CR1-AR1*CX ) 

V=6X*CX*< AR1*DR2-AR2*DR1)+BR1*DX*( AR2*CX-AX*CR2) 

P=AX**2*(CR1*DR2-CR2*DR1)+AX*AR2*(CX*DR1-CR1*DX)+AR2*BX*CX*CR1 

Q=CX**2*( AR1*BR2-AR2*BR1)+CX*CR2*( AX*BR1-AR 1*BX ) -AX*BR2*CX*CR 1 

C=(A+B)/(D+E) 

HX=AM1*C**3+BM1*C**2+CM1*C+DM1 

G  X= ( A+B ) / ( D+  E ) - ( U  + V ) / ( P+0 ) 

WX=C**2*( AM1R1*CX-AX*CM1R1+AM1*CR1+AR1*CM1 ) +C* ( AM1*DR 1+AR1*DR 1+BR 1 
1*CM1+BR1*CR1+AM1R1*DX-AX*DM1R1+BM1R1*CX-BX*CM1R1 ) +BM1*DR 1+BR1*DM1 + 
2BM1R1*DX-BX*DM1R1 

WY=2.*(CX*C+DX)*IAM1*C+BM1 ) * ( AR 1*C+BR1 ) / ( AX*C+BX ) 

WZ=2.*C*( AX*C+BX)**2-AX*DX+BX*CX 

C12=(WX-WY)/WZ 

12  PRINT  2.C.XM.R 
PRINT  4,R,GX,HX 

13  PUNCH  3.C.XM.C12 
PUNCH  4.R.GX.HX 
GO  TO  10 

END 
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INPUT  DATA  CF  PRC 

!GRAM  1. 

5 

.00200E-00 
+.001     +.001 

0.50200E-00 

5 

.01200E-00 
■1.01      +.01 

0.50000E-00 

5 

•00100E-00 
-  .0001   -0.0001 

0.50200E-00 

5 

.35400E-00 
+.001     +.001 

O.OAOOOE-OO 

5 

.37000E-00 
.005     0.01 

0.02100E-00 

5 

.32450E-00 
-  .0001   +0.0 

O.17900E-00 

5 

•32450E-00 
+  .0      +0.01 

0.17900E-00 

1 

.5O2OOE-0O 
+  .0      +0.U 

0.00100E-00 

1 

.34390E-00 
+  .0      +0.0 

0.03790E-00 

3 

•35450E-00 
-  .005    +0.001 

0.04250E-00 

3 

•33450E-C0 

-  .005    -0.^05 

0.02790E-00 

3 

•32420E-00 
+  .01     +0. 

0.01790E-00 
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COMPUTER  RESULTS  CF 

PRCGRAM  1. 

C" 

5614.91378E+11 

XM  =  2.0UOOOE-03 

C12=-3.87070E-06 

R  = 

5u20.00uOGE-04 

GX  =  7.88194E-09 

HX=  2.98650E-10 

C  = 

Z0C1.26781E-10 

XM=  3.00000E-03 

C12=-2.06148E-05 

R  = 

503C.0OUO0E-04 

GX  =  2.98780E-08 

HX=  1.66781E-09 

C" 

4998.09563E-10 

XM=  4.00000E-03 

C12=-6.02374E-05 

R  = 

5040.00uOOE-04 

GX=  7.72587E-08 

HX=  5.77638E-09 

C  = 

1026.14803E-U9 

XM=  5.00000E-03 

C12=-1.34925E-04 

R" 

5u50.00UOOE-04 

GX  =  1.60096E-07 

HX=  1.53600E-08 

C  = 

1859.60756E-09 

XM=  6.0G000E-03 

C12=-2.58965E-04 

R  = 

5060.0O000E-04 

GX=  2.85O99E-07 

HX=  3.44926E-08 

C  = 

1946.70853E-08 

XM=  1.2000CE-02 

C12=-3.34778E-03 

R  = 

500O.00O00E-G4 

GX=  8.25184E-06 

HX=  7.72637E-07 

C« 

1462.97428E-07 

XM=  2.2UOOOE-02 

C12=-3.19896E-02 

R  = 

51CG.00u0uE-04 

GX=-3.42133E-05 

HX=  1.10602E-05 

C  = 

4810.B4U67E-G7 

XM=  3.200G0E-02 

C12=-1.58675E-01 

R  = 

5200.00000E-G4 

GX=  8.93279E-04 

HX=  4.20783E-05 

C" 

1112.07029E-06 

XM=  4.2000CE-02 

C12=-1.25363E-00 

R  = 

5300.00GOOE-04 

GX=  1.38002E-03 

HX=  3.95923E-05 

C  = 

2173.34047E-06 

XM=  5.20000E-02 

C12  =  9.15740E-01 

R  = 

54G0.00G00E-04 

GX=  2.43379E-03 

HX=-3.32558E-04 

C  = 

66C1.89946E-12 

XM=  1.000UOE-03 

C12=  2.03680E-07 

R  = 

5020.00G00E-04 

GX=  6.41334E-1U 

HX=  1.68023E-1! 

C  = 

4784.87552E+12 

XM=  9.00000E-04 

C12=  2.53433E-07 

R  = 

5019.00000E-04 

GX=  4.58203E-10 

HX=  1.09037E-11 

5k 

c= 

1054.58622E-05 

XM  = 

3.54000E-01 

C12=-1.02713E+01 

R  = 

4000.00000E-05 

GX  = 

'  4.49330E-02 

HX=-2.80293E-01 

C  = 

1062.72215E-U5 

XM  = 

3.550OOE-01 

C12=-9.23096E-00 

R  = 

410G.GOGGOE-U5 

GX  = 

i  4.46904E-02 

HX=-2.89962E-01 

C  = 

1G70.46877E-05 

XM  = 

3.560G0E-O1 

C12=-8.32577E-00 

R  = 

4200.00000E-05 

GX  = 

:  4.44540E-02 

HX=-2.99867E-01 

C« 

1077.83685E-05 

XM  = 

3.57000E-01 

C12=-7.53435E-00 

R  = 

43GG.00GG0E-G5 

GX  = 

■  4.42231E-02 

HX=-3.10012E-01 

C" 

1G84.83631E-05 

XM  = 

3.58000E-01 

C12=-6.83936E-00 

R  = 

440G.00GOGE-G5 

GX  = 

i  4.39971E-02 

HX=-3.20401E-01 

C  = 

6045.59426E-06 

XM  = 

3.700GOE-01 

C12=-1.50812E+01 

R  = 

2100.00000E-05 

GX  = 

1  4.49355E-02 

HX=-2.40059E-01 

C" 

7566.30695E-G6 

XM  = 

3.750G0E-01 

C12=-6.99655E-00 

R« 

3100.00UOUE-05 

GX  = 

=  4.12903E-02 

HX=-3.18245E-01 

C  = 

8743.54689E-U6 

XM  = 

3.8G0GGE-01 

C12=-3.72768E-00 

R  = 

410G.G0UOGE-U5 

GX  = 

>  3.96779E-02 

HX=-4.11281E-01 

C  = 

9642.64906E-06 

XM  = 

3.8500GE-01 

C12=-2.19897E-00 

R  = 

51O0.0OO0OE-05 

GX  = 

=  3.86665E-02 

HX=-5.20351E-01 

C  = 

1030.70943E-05 

XM  = 

3.900G0E-01 

C12=-1.42012E-00 

R« 

6100.00OG0E-05 

GX  = 

>  3.78118E-02 

HX=-6.46227E-01 

C  = 

1054.58622E-05 

XM  = 

3.540GOE-01 

C12=-1.02713E+01 

R  = 

40OG.0OG0OE-G5 

GX  = 

>  4.49330E-02 

HX=-2.80293E-01 

C" 

1053.75083E-05 

XM  = 

3.53900E-01 

C12=-1.03838E+01 

R« 

3990.00000E-G5 

GX  = 

=  4.49577E-02 

HX=-2.79339E-01 

C" 

1075.54104E-05 

XM  = 

3.48900E-01 

C12=-1.42297E+01 

R  = 

3890.00G00E-05 

GX  = 

=  4.62912E-02 

HX=-2.53947E-01 
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c= 

7864.13244E-06 

XM  = 

3.24500E-01 

C12=  9.33546E+01 

R  = 

1790.00000E-C5 

GX  = 

■  5.84102E-02 

HX=-1.05281E-01 

C» 

8157.52417E-06 

XM  = 

3.24500E-01 

C12  =  9.58772E+01 

R  = 

189C.00GOOE-05 

GX  = 

•  5.78839E-02 

HX=-1.08110E-01 

C" 

8441.94571E-06 

XM  = 

3.24500E-01 

C12=  9.89806E+01 

R* 

1990.00UOOE-U5 

GX  = 

i  5.74477E-02 

HX=-1.10957E-01 

C" 

8717.59689E-06 

XM  = 

3.24500E-01 

C12"  1.02788E+02 

R« 

2090.00000E-05 

GX  = 

■  5.70875E-02 

HX="1.13822E-01 

C" 

8984.67593E-06 

XM  = 

3.24500E-01 

C12=  1.07471E+02 

R  = 

2190.00G00E-G5 

GX  = 

>  5.67914E-02 

HX=-1.16706E-01 

C  = 

2044.28003E-06 

XM  = 

3.24500E-00 

C12=-3.81504E-03 

R* 

179C.00000E-05 

GX  = 

>  1.25790E-02 

HX=-6.70091E+05 

C  = 

2044.29912E-06 

XM  = 

3.24400E-00 

C12=-3.81791E-03 

R* 

1790.00000E-05 

GX  = 

■  1.25817E-02 

HX=-6.68566E+05 

C" 

2044.31909E-06 

XM  = 

3.243U0E-00 

C12=-3.82078E-03 

R« 

1790.00000E-05 

GX  = 

'  1.25845E-02 

HX=-6.67045E+05 

C  = 

2044.33876E-06 

XM  = 

3.24200E-00 

C12=-3.82366E-03 

R  = 

1790.00000E-O5 

GX  = 

'  1.25873E-02 

HX=-6.65526E+05 

C  = 

2044.35791E-C6 

XM  = 

3.24100E-00 

C12=-3.82653E-03 

R  = 

1790.00000E-05 

GX  = 

■  1.25901E-02 

HX=-6.64011E+05 

C  = 

1977.08684E-06 

XM  = 

5.02000E-01 

C12=-1.29676E-00 

R  = 

IUOU.OO0GUE-O6 

GX  = 

=-9.60148E-02 

HX=-1.18903E-00 
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c= 

1096.18200E-05 

XM  = 

3.43900E-01 

C12=-2.04427E+01 

Ra 

379Q.00OO0E-05 

GX  = 

■  4.77748E-02 

HX=-2.30236E-01 

c= 

1089.69948E-05 

XM  = 

3.54500E-01 

C12=-8.62141E-00 

R  = 

4250.00000E-05 

GX  = 

■  4.48423E-02 

HX=-2.96020E-01 

c» 

1142.95223E-05 

XM  = 

3.49500E-01 

C12=-1.01324E+01 

Re 

4350.00UOOE-05 

GX  = 

>  4.62297E-02 

HX=-2.79586E-01 

C  = 

1194.26414E-05 

XM  = 

3.44500E-01 

C12=-1.20154E+01 

R  = 

4450.00000E-05 

GX  = 

>  4.77956E-02 

HX=-2.63286E-01 

C  = 

1091.75240E-05 

XM  = 

3.44500E-01 

C12=-1.96992E+01 

Ra 

3790.00000E-05 

GX  = 

■■    4.75856E-02 

HX=-2.32447E-01 

C  = 

1038.67725E-05 

XM  = 

3.39500E-01 

C12=-4.54421E+01 

Ra 

3290.00000E-05 

GX  = 

•  4.92179E-02 

HX=-1.93753E-01 

C  = 

9726.00035E-06 

XM  = 

3.34500E-01 

C12=-2.29778E+02 

Ra 

279G.00000E-05 

GX  = 

■■    5.12844E-02 

HX=-1.59976E-01 

C  = 

8902.15514E-06 

XM  = 

3.29500E-01 

C12=  1.75740E+02 

R  = 

2290.00000E-C5 

GX  = 

■■    5.41146E-02 

HX=-1.30648E-01 

O 

7864.13244E-06 

XM  = 

3.24500E-01 

C12=  9.33546E+01 

Ra 

1790.00000E-05 

GX  = 

•  5.841C2E-02 

HX=-1.05281E-01 

C  = 

7883.41285E-C6 

XM  = 

3.24200E-01 

C12=  9.24051E+01 

Ra 

1790.00000E-05 

GX  = 

■■    5.85658E-02 

HX=-1.04715E-01 

C  = 

8535.95669E-06 

XM  = 

3.14200E-01 

C12  =  7.45334E+01 

Ra 

179C.OOOOOE-05 

GX  = 

■  6.51697E-02 

HX=-8.70843E-02 

C  = 

9185.78979E-06 

XM  = 

3.04200E-01 

C12=  6.74310E+01 

R  = 

1790.00000E-05 

GX  = 

'  7.63863E-02 

HX=-7.16546E-02 
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PROGRAM  2.  CALCULATION  OF  THE  SPECTRUM  OF  AN  EQUAL  AMPLITUDES 
NOTCH  FILTER  OF  THE  GOLDMAN  TYPE 


C   EQUAL  AMPLITUDES  NOTCH  FILTER  OF  THE  GOLDMAN  TYPE 

1  READ,  wl,  W2,  DELW 
PUNCH  4 

4  FORMATdlX,  5H0MEGA1UX.4HEVAL/1 
3  W  =  Wl 

X  =  ( 1.0-0.41421*W**2)**2+(0.29  2895*W-0.121315*W**3)**2 

Y=(1.0-1.82842*W**2)**2+(3.707105*W-C.121315*W**3)**2 

EVAL=X/Y 

PUNCH  b     ,  W   ,  EVAL 

5  FORMAT     (     5X,    F12.6     ,    5X.    F12.6) 
Wl    =    Wl    +    DELW 

IF     (     W2    -    W     )     2.3.3 

2  GO    TO    1 
END 
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C   C   DATA  FCR  GRAPH  OF 

VALUE  VS  OMEGA  -  EQUAL 

AMPLITUDES 

NOTCH  FILTER 

OF  THt  GOLDMAN  TYPt 

OMEGA 

EVAL 

.000000 

1.0000  0  0 

.050000 

.973581 

.100000 

.901446 

.150000 

.800658 

.2000CC 

.689567 

.250000 

.581728 

.300000 

.484553 

.350000 

.400739 

.4J0ut.u 

.330231 

.450000 

.271696 

.5U0OOO 

.223398 

.550000 

.183617 

.60OOOO 

.150833 

.650000 

.123766 

.700000 

.101366 

.750000 

.082787 

.800000 

.067348 

.850000 

.054500 

.900000 

.043804 

.950000 

.034903 

1*000000 

.027509 

1.050000 

.021386 

1.100000 

.016340 

1.150000 

.012211 

1.200000 

.008867 

1.250000 

.006196 

1.300000 

.004106 

1.350000 

.002518 

1.400000 

.001366 

1.450000 

.000594 

1.500000 

.000152 

1.550000 

.000000 

1.600000 

.000103 

1.650000 

.000429 

1. 700000 

.000952 

1.750000 

.001648 

1.800000 

.002499 

1.850000 

.003485 

1.900000 

.004593 

1.950000 

.005809 

2.000000 

.007120 

2.050000 

.008517 

2.100000 

.009992 

2.15oOuo 

.011535 

2.2oo0ou 

.013140 

2.250000 

.014801 

2.300000 

.016513 

2.350000 

.018270 

2.400000 

.020069 

2.450000 

.021905 
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OMEGA 

EVAL 

2.  5^0000 

.023777 

2.550000 

.025679 

2.600000 

.027611 

2.650000 

.029570 

2.700000 

.031554 

2.750000 

.033560 

2.800000 

.035588 

2.850000 

.037637 

2.900000 

.039704 

2.950000 

.041788 

3.0u000u 

.043890 

3.050000 

.046007 

3.100000 

.048139 

3.150000 

.050286 

3.200000 

.052447 

3.250000 

.054621 

3.300000 

.056808 

3.350000 

.059007 

3.400000 

.061217 

3.450000 

.063440 

3.500000 

.065674 

3.550000 

.067918 

3.6u0000 

.070173 

3.65UUUJ 

.072439 

3.700000 

.074715 

3.750000 

.077000 

3.800000 

.079296 

3.850000 

.081601 

3.900000 

.083915 

3.950000 

.086239 

4.000000 

.088572 

4.050000 

.090913 

4. 100000 

.093263 

4.150000 

.095622 

4.200000 

.097989 

4.250000 

.100365 

4.3^0000 

.102748 

4.350000 

.105140 

4.400000 

.107539 

4.450000 

.109946 

4.500000 

.112361 

4.550000 

.114783 

4.600000 

.117213 

4.650000 

.119649 

4.700000 

.122093 

4.750000 

.124543 

4.800000 

.127000 

4.850000 

.129464 

4.900000 

.131934 

4.950000 

.134410 

5.OUO0OO 

.136893 
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0ME6A 

FVAL 

.000000 

1.000000 

1.000000 

.027509 

2.000000 

.007120 

3.000000 

.043890 

4. OoOOOO 

.088572 

5.000000 

.136893 

6. OoOOOO 

.187582 

7.000000 

.239365 

8.0O0UOU 

.290978 

9.000000 

.341341 

10.000000 

.389635 

11.000000 

.435303 

12.000000 

.478020 

13.000000 

.517647 

14.000000 

.554178 

15.000000 

.587703 

16.000000 

.618370 

17.000000 

.646365 

18. OoOOOO 

.671889 

19.000000 

.695146 

20.000000 

.716338 

21.000000 

.735654 

22.000000 

.753274 

23.000000 

.769362 

24.000000 

.784066 

25.000000 

.797525 

26.000000 

.809859 

27.000000 

.821179 

28.000000 

.831584 

29.0oOooO 

.841162 

30.000000 

.849992 

31.000000 

.858146 

32.000000 

.865686 

33.000000 

.872668 

34.000000 

.879145 

35.000000 

.885160 

36.000000 

.890755 

37.000000 

.895966 

38.000000 

.900826 

39. OOOOOO 

.905366 

40.0o0000 

.909610 

41.000000 

.913584 

42.000000 

.917310 

43.000000 

.920806 

44.000000 

.924091 

45.000000 

.927181 

46.000000 

.930090 

47.000000 

.932833 

48.000000 

.935422 

49.000000 

.937866 

50.000000 

.940178 
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CMEGA 

EVAL 

50.0uuoou 

.940178 

60.000000 

.957672 

70.0U00Ou 

.968544 

80.000000 

.975735 

90.000000 

.980728 

100.000000 

.984332 

] 10.000000 

.987015 

120.000000 

.989066 

130.000000 

.990668 

140.000000 

.991943 

150.0U0uGu 

.992974 

160.000000 

.993820 

170.0u0000 

.994522 

180.0u0000 

.995110 

190.000000 

.995609 

200.000000 

.996036 

210.000000 

.996403 

220.000000 

.996721 

230.000000 

.996999 

240.000000 

.997244 

250.000000 

.997459 

260.000000 

.997650 

270.000000 

.997821 

280.000000 

.997973 

290.000000 

.998110 

300.000000 

.998234 

310.000000 

.998346 

320.000000 

.998448 

330.000000 

.998540 

340.000000 

.998625 

350.000000 

.998702 

360.000000 

.998773 

370.000000 

.998838 

380.000000 

.998899 

390.000000 

.998954 

400.0UOO00 

.999006 

410»QOO0O0 

.999054 

420.000000 

.999098 

430.000000 

.999140 

440.000000 

.999179 

450.000000 

.999215 

460.000000 

.999248 

470.000000 

.999280 

480.000000 

.999310 

490.000000 

.999337 

500.000000 

.999364 
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OMEGA 

EVAL 

500.000000 

.999364 

600.000000 

.999558 

700. 000000 

.999675 

800.000000 

.999751 

900.000000 

.999803 

lOOO.OOuOGO 

.999841 

1100.0i>0000 

.999868 

1200.000000 

.999889 
.999906 

1300.000000 

1400.000000 

.999919 

1500.000000 

.999929 

1600.000000 

.999938 

1700.000000 

.999945 

1800.000000 

.999951 

1900.000000 

.999956 

2000.000000 

.999960 

210O.OuOOO0 

.999964 

2200.0o000G 

.999967 

2300.000000 

.999970 

2400.000000 

.999972 

2500.000000 

.999975 

2600.000000 

.999976 

2700.000000 

.999978 

2800.000000 

.999980 

2900.000000 

.999981 

3000.000000 

.999982 

3100.000000 

.999983 

3200.000000 

.999984 

3300.000000 

.999985 

3400.0O00G0 

.999986 

3300.000000 

.999987 

36OG.0U00OO 

.999988 

3700.000000 

.999988 

3800.000000 

.999989 

3900.000000 

.999990 

4000.000000 

.999990 

4100.000000 

.999991 

4200.000000 

.999991 

4300.000000 

.999991 

4400.000000 

.999992 

4500.000000 

.999992 

4600.000000 

.999993 

4700.000000 

.999993 

480O.0OUO00 

.999993 

4900.000000 

.999993 

5000.000000 

.999994 

5100.000000 

.999994 

5200.000000 

.999994 

5300.000000 

.999994 

5400.000000 

.999995 
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OMEGA 

EVAL 

5500*000000 

.999995 

5600.000000 

.999995 

5700.000000 

.999995 

5800.0^0000 

.999995 

5900.000000 

.999995 

6000.000000 

.999996 

6100.000000 

.999996 

6200.000000 

.999996 

6300.000000 

.999996 

6400.000000 

.999996 

6500.000000 

.999996 

6600.000000 

.999996 

6700.000000 

.999996 

6800.000000 

.999997 

6900.000000 

.999997 

7000.000000 

.999997 

7100.000000 

.999997 

7200.000000 

.999997 

7300.000000 

.999997 

7'*O0. 000000 

.999997 

7500.000000 

.999997 

7500.000000 

.999997 

7700.000000 

.999997 

7800.000000 

.999997 

7900.000000 

.999997 

8000.000000 

.999998 

8100.000000 

.999998 

8200.0O0000 

.999998 

8300.000000 

.999998 

8400.000000 

.999998 

8500.000000 

.999998 

8600.000000 

.999998 

8700.000000 

.999998 

8800.000000 

.999998 

8900.000000 

.999998 

9000.000000 

.999998 

9100.000000 

.999998 

9200.000000 

.999998 

9300.000000 

.999998 

9400.000000 

.999998 

9500.000000 

.999998 

9600.000000 

.999998 

97C0. 000000 

.999998 

9800.000000 

.999998 

9900.000000 

.999998 

10000.000000 

.999998 

10100.000000 

.999998 

10200.000000 

.999998 

10300.000000 

.999998 

10400.000000 

.999999 

10500.000000 

.999999 

^ 


PRCGRAM  3.  CALCULATION  OF  THE  SPECTRUM  CF  AN  UNEQUAL 
AMPLITUDES  NOTCH  FILTER 


C   UNEQUAL  AMPLITUDES  NOTCH  FILTER 

1  READ.  Wl,  W2.  DELW 
PUNCH  4 

A  FORMATU1X.  5HOMEGA1 CX .4HEVAL/ ) 
3  W  =  Wl 

X  =  144.*W**6  +  2  766.25*W**<t  +  5172.*W**2  +  256. 

EVAL=(2.25*W**6-2.*W**4-92.*W**2+2  56. )/X 

PUNCH  5  .  W   .  EVAL 
5  FORMAT  (  5X,  F12.6  ,  5X  .  F12.6) 

Wl  =  Wl  +  DELW 

IF  (  W2  -  W  )  2.3.3 

2  GO  TO  1 
END 
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C   C   DATA  FCR  GRAPH 

OF  VALUE  VS  CMEGA  -  UNEQUAL 

AMPLITUDES 

NOTCH  FILTEk 

CMEGA 

EVAL 

.000000 

1.000000 

.050000 

.951004 

.100000 

.828190 

.150000 

.679369 

.200000 

.539929 

.250000 

.424075 

.300000 

.332941 

.350000 

.262743 

.4^0000 

.208878 

.450000 

.167358 

.500000 

.135088 

.550000 

.109763 

.600000 

.089695 

.650000 

.073644 

.700000 

.060697 

.750000 

.050175 

.800000 

.041565 

.850000 

.034480 

.900000 

.028621 

.950000 

.023754 

1.000000 

.019698 

1.050000 

.016309 

1.100000 

.013470 

1.150000 

.011089 

1.200000 

.009091 

1.250000 

.007414 

1.300000 

.006006 

1.350000 

.004828 

1.400000 

.003843 

1.450000 

.003022 

1.500000 

.002342 

1.550000 

.001783 

1.600000 

.001326 

1.650000 

.000957 

1.7^0000 

.000664 

1.750000 

.000436 

1.800000 

.000264 

1.850000 

.000141 

1.900000 

.000059 

1.950000 

.000014 

2.000000 

0.000000 

2.050000 

.000013 

2.100000 

.000049 

2.150000 

.000105 

2.2OO000 

.000178 

2.250000 

.000267 

2.300000 

.000368 

2.350000 

.000479 

2.400000 

.000601 

2.450000 

.000730 
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OMEGA 

EVAL 

2.500000 

.000866 

2.550000 

.001007 

2.600000 

.001153 

2.650000 

.001304 

2.700000 

.001457 

2.750000 

.001612 

2.800000 

.001770 

2.850000 

.001929 

2.9^0000 

.002089 

2.95000^, 

.002250 

3.000000 

.002411 

3.050000 

.002572 

3.100000 

.002733 

3.150000 

.002893 

3.200000 

.003053 

3.250000 

.003212 

3.300000 

.003370 

3.3500^u 

.003526 

3.4OG00U 

.003682 

3.450000 

.003836 

3.500000 

.003988 

3.550000 

.004139 

3.6o0000 

.004289 

3.650000 

.004437 

3.700000 

.004583 

3.750000 

.004727 

3.800000 

.004870 

3.850000 

.00501C 

3.900000 

.005150 

3.950000 

.005287 

4.0uOOOO 

•  0C542i: 

4.050000 

.005556 

4.100000 

.005688 

4.150000 

.005817 

4.200000 

.005945 

4.250000 

.006072 

4.300000 

.006196 

4.350000 

.006319 

4.400000 

.006440 

4.450000 

.006559 

4.500000 

.006676 

4.55000C 

.006791 

4.6U0000 

.006905 

4.650000 

.007017 

4.700000 

.007128 

4.750000 

.007236 

4.800000 

.007343 

4.850000 

.007449 

4.900000 

.007553 

4.950000 

.007655 

5.000000 

.007755 
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OMEGA 

EVAL 

.000000 

1.000000 

1.000000 

.019698 

2.000000 

0.000000 

3.000000 

.002411 

4.000000 

.005422 

5.000000 

.007755 

6.000000 

.009469 

7.000000 

.010725 

8.O0OOOO 

.011657 

9.0o0000 

.012360 

lU.OoOOUO 

.012900 

11.000000 

.013320 

12.000000 

.013653 

13.000000 

.013921 

14.000000 

.014139 

15.000000 

.014318 

16.000000 

.014468 

17.000000 

.014593 

18.000000 

.014700 

19.000000 

.014791 

2u.0u0000 

.014869 

21.000000 

.014937 

22.000000 

.014996 

23.000000 

.015048 

24.000000 

.015094 

25.000000 

.015135 

26.000000 

.015171 

27.000000 

.015203 

28.000000 

.015232 

29.000000 

.015258 

30.000000 

.015282 

31.000000 

.015303 

32.000000 

.015323 

33.000000 

.015341 

34.000000 

.015357 

35.000000 

.015372 

36.000000 

.015386 

37.000000 

.015398 

38.000000 

.015410 

39.000000 

.015421 

40.000000 

.015431 

41.000000 

.015440 

42.000000 

.015449 

43.000000 

.015457 

44.0OC000 

.015464 

45.000000 

.015471 

46.000000 

.015478 

47.0o0000 

.015484 

48.000000 

.015490 

49.000000 

.015495 

50.000000 

.015500 
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OMEGA 

EVAU 

50.000000 

.015500 

60.000000 

.015538 

70.000000 

.015561 

80.000000 

.015576 

90.000000 

.015586 

100.000000 

.015594 

110.000000 

.015599 

120.000000 

.015603 

130.000000 

.015606 

140.000000 

.015609 

150.000000 

.015611 

160.000000 

.015613 

170.000000 

.015614 

180.000000 

.015615 

19C. 000000 

.015616 

200.000000 

.015617 

210.000000 

.015618 

220.000000 

.015619 

230.000000 

.015619 

240.000000 

.015620 

250.000000 

.015620 

260.000000 

.015620 

270.000000 

.015621 

280.000000 

.015621 

290.0O0000 

.015621 

300.000000 

.015622 

310.000000 

.015622 

320.000000 

.015622 

330.000000 

.015622 

340.000000 

.015622 

350.000000 

.015622 

360.000000 

.015623 

370.000000 

.015623 

380.000000 

.015623 

390.000000 

.015623 

400.000000 

.015623 

410.000000 

.015623 

420.000000 

.015623 

430.000000 

.015623 

440.000000 

.015623 

450.0U0000 

.015623 

460.000000 

.015624 

470.000000 

.015624 

480.000000 

.015624 

490.000000 

.015624 

500.000000 

.015624 

510.000000 

.015624 
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CMFGA 

FVAL 

500.000000 

.015624 

600.000000 

.015624 

700.000000 

.015624 

BOO.OuOuOC 

.015625 

900.0^0000 

.015625 

lOOO.OuOOOO 

.015625 

llOO.OOuOOU 

.015625 

1200.000000 

.015625 

1300.000000 

.015625 

1400.000000 

.015625 

1500.000000 

.015625 

1600.000000 

.015625 

1700.000000 

.015625 

1800.000000 

.015625 

1900.000000 

.015625 

2000.000000 

.015625 

2100.000000 

.015625 

2200.0uOOOO 

.015625 

2300.000000 

.015625 

2400.000000 

.015625 

2500.000000 

.015625 

2600.000000 

.015625 

2700.000000 

.015625 

2800.000000 

.015625 

2900.000000 

.015625 

3000.000000 

.015625 

3100.000000 

.015625 

3200.0uuGUO 

.015625 

330G.0u0uuu 

.015625 

3400.0UOUCO 

.015625 

3500.000000 

.015625 

3600.0U0000 

.015625 

3700.000000 

.015625 

3800.000000 

.015625 

3900.000000 

.015625 

4000.000000 

.015625 

4100.000000 

.015625 

4200.000000 

.015625 

4300.0^0000 

.015625 

4400.000000 

.015625 

4500.000000 

.015625 

4600.000000 

.015625 

4700.000000 

.015625 

4800.000000 

.015625 

4900.000000 

.015625 

5000.000000 

.015625 

5100.000000 

.015625 

5200.000000 

.015625 

5300.000000 

.015625 

5400.000000 

.015625 
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PROGRAM  4.  CALCULATION  OF  THE  SPECTRUM  OF  AN  EQUAL 
AMPLITUDES  NOTCH  FILTER 


C   EQUAL  AMPLITUDES  NOTCH  FILTER 

1  READ9.W1.  W2.  DELW.  A 
9  F0RMATI3E10.4.F2.0) 

IF  (  A  )  3.3.6 

3  W  =  Wl 

4  F0RMATI11X.  5H0MEGA10X .4HE VAL/ I 

Y  =  0.5999928*W**4-28.7633238*W**2+322. 3679696 
Z  =  0.03053742*W**5-3.8378482*W**3+165.084187*W 
X  ■  7.4503121»W**4-677.1501976*W**2+322. 3679696 
P  =  G.u3053742*W**5-122.384064*W**3+1291.817275*W 
EVAL  =  (Y**2+Z»*2)/(X**2+P**2 ) 
PUNCH  5  ,  W   »  EVAL 
J  FORMAT  (  5X.  F12.6  .  5X.  E14.6  ) 
Wl  =  Wl  +  DELW 
IF  (  W2  -  0  )  2.3.3 

2  GO  TO  1 
6  STOP 

END 
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C   C   DATA  FOR  GRAPH  : 

3F  VALUE  VS  CMEGA  -  EQUAL 

AMPLILUOES 

NOTCH  FILTER 

CMEGA 

EVAL 

.OOOOOO 

0.100000E+C1 

.050000 

0.971406 

.100000 

0.894633 

.150000 

0. 790440 

,2oOOOo 

0. 679522 

.250000 

0.575533 

i3OOOO0 

0.484708 

.350000 

0.408377 

.400000 

0.345445 

.450000 

0.293955 

.500000 

0.251858 

.550000 

0.217331 

.600000 

0.188859 

.650000 

0.165226 

.7U0000 

0.145470 

.750000 

0.128837 

•SoUOOO 

0.114735 

.850000 

0.102696 

•9o0000 

0.923526E-01 

.950000 

0.834105E-01 

l.oooooo 

0.756351E-01 

1.050000 

0.688373E-01 

1.100000 

0.628638E-01 

1.150000 

0.575894E-01 

1.200000 

0.529111E-01 

1.250000 

0.487442E-01 

1.300000 

0.450179E-01 

1.3500C0 

0.416732E-C1 

1.400000 

0.386606E-01 

1.450000 

0.359380E-01 

1.500000 

0.334699E-01 

1.550000 

O.312259E-01 

1.600000 

0.291800E-01 

1.650000 

0.273098E-01 

1.700000 

0.255960E-01 

1.750000 

0.240219E-01 

1.800000 

0.225727E-01 

1.850000 

0.212359E-01 

1.900000 

0. 2000016-01 

1.950000 

0.188557E-01 

2.0o0000 

0.177938E-01 

2.050000 

0.168069E-01 

2.1o0000 

0.158881E-01 

2.150000 

0.150315E-01 

2.200000 

0.142316E-01 

2.250000 

0.134836E-01 

2.300000 

0.127832E-U1 

2.350000 

0.121265E-C1 

2.400000 

0.115100E-01 

2.450000 

0.109307E-01 
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2.500000 

0.103855E-01 

2.550000 

0.987206E-02 

2.600000 

0.938794E-02 

2.650000 

0.893103E-02 

2.7^0000 

0.849941E-02 

2.750000 

0.809131E-02 

2.800000 

0.770513E-02 

2.850000 

0.733939E-02 

2.9^0000 

0.699275E-02 

2.950000 

0.666398E-02 

3.000000 

0.635192E-02 

3.050000 

0.605555E-02 

3.100000 

0.577389E-02 

3.150000 

0.550605E-C2 

3.200000 

0.525122E-02 

3.250000 

0.500863E-02 

3.300000 

0.477758E-02 

3.350000 

0.455742E-02 

3.400000 

0.434754E-02 

3.450000 

0.414737E-02 

3.500000 

0.395640E-02 

3.550000 

0.377413E-02 

3.600000 

0.360010E-02 

3.650000 

0.343391E-02 

3.700000 

0.327514E-02 

3.750000 

0.312343E-02 

3.800000 

0.297842E-02 

3.850000 

0.283980E-C2 

3.900000 

0.270725E-02 

3.950000 

0.258049E-02 

4.000000 

0.245926E-02 

4.050000 

0.234329E-02 

4.1u0000 

0.223234E-02 

4.150000 

0.212621E-02 

4.200000 

0.202466E-02 

4.250000 

0.192750E-02 

4.300000 

0.183455E-02 

4.350000 

0.174562E-02 

4.400000 

0.166054E-02 

4.450000 

0.157916E-02 

4.500000 

0.150133E-02 

4.550000 

0.142689E-02 

4.600000 

0.135572E-02 

4.650000 

0.128769E-02 

4.700000 

0.122266E-02 

4.750000 

0.116054E-02 

4.800000 

0.110120E-C2 

4.850000 

0.104454E-02 

4.900000 

0.990459E-03 

4.950000 

0.938864E-03 

5.000000 

0.889661E-L3 
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5.05000C 

0.842765E-03 

5.100000 

0.798091E-03 

5.150000 

0.755561E-03 

5.200000 

0.715099E-C3 

5.250JOO 

0.676633E-03 

5.300000 

0.640094E-C3 

5.350000 

0.605414E-03 

5.400000 

0.572531E-03 

5.450000 

0.541385E-03 

5.500000 

0.5U916E-03 

5.550000 

0.48407QE-03 

5.600000 

0.457792E-03 

5.650000 

0.433031E-03 

5.700000 

0.409739E-03 

5.750000 

0.387867E-03 

5.800000 

0.367370E-03 

5.850000 

0.348206E-03 

5.900000 

0.330331E-03 

5.9500OL 

0.313705E-03 

6.000000 

0.298290E-03 

6.050000 

0.284048E-03 

6.100000 

0.270943E-03 

6.150000 

0.258941E-03 

6.200000 

0.248008E-03 

6.250000 

0.238111E-03 

6.300000 

0.229221E-03 

6.350000 

0.221307E-03 

6.400000 

0.214341E-03 

6.45000Q 

0.208294E-03 

6.500000 

0.203139E-03 

6.550000 

0.198852E-03 

6.6UU000 

0.195407E-03 

6.650000 

0.192779E-03 

6.700000 

0.190946E-03 

6.750000 

0.189885E-03 

6.800000 

0.189574E-03 

6.850000 

0.189992E-03 

6.900000 

0.191120E-03 

6.950000 

0.192936E-03 

7.000000 

0.195423E-03 

7.050000 

0.198562E-03 

7.10000C 

0.202336E-03 

7.150000 

0.206726E-03 

7.2^0000 

0.211717E-03 

7.250000 

0.217292E-03 

7.300000 

0.223436E-03 

7.35000C 

0.230133E-03 

7.400000 

0.237370E-03 

7.450000 

0.245131E-03 

7.500000 

0.253404E-03 

7.550000 

0.262175E-03 

7.600000 

0.271431E-03 
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7.650000  0.281160E-03 

7.700000  0.291349E-03 

7.750000  0.301988E-U3 

7.80000U  0.313064E-03 

7.850000  0.324567E-03 

7.900000  0.336487E-C3 

7.950000  0.348812E-03 

B. 000000  0.361533E-Q3 

8.05uuuo  0.374639E-03 

8.1uu0uu  0.388123E-C3 

8.150000  0.401975E-03 

8.200000  0.416185E-03 

8.250000  0.430745E-03 

8.300000  0.445648E-C3 

8.350000  0.460884E-03 

8.400000  0.476446E-03 

8.450000  0.492327E-03 

8.50000C  0.508519E-03 

8.55uuou  0.525015E-03 

8.600000  0.541808E-03 

8.650000  0.558892E-03 

8.7O0000  0.576259E-03 

8.750UOO  0.593904E-03 

8.8O0OU0  0.611820E-03 

8.850000  0.630002E-03 

8.900000  0.648443E-03 

8.950000  0.667138E-C3 

9.000000  0.686082E-C3 

9.050000  0.70527UE-Q3 

9.100000  0.724695E-03 

9.150000  0.744354E-03 

9.2U0000  0.764241E-03 

9.250000  0.784352E-03 

9. 300000  0.804681E-03 

9«350000  0.825226E-03 

9.4uu000  0.845981E-03 

9.450000  0.866942E-03 

9.5^0000  0.888106E-03 

9.550000  0.909468E-C3 

9.60000C  0.931024E-03 

9.650000  0.952771E-03 

9.700000  0.9747C6E-C3 

9.750000  0.996824E-03 

9.8OU0OO  0.101912E-02 

9.850O0O  0.104160E-02 

9.9^0000  0.106425E-02 

9.950000  0.108707E-02 

lO.OOOOuO  0.111006E-02 

10.050000  0.113321E-02 
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10.000000 

0.111006E-02 

11.000000 

0.160113E-C2 

12.000000 

0.214152E-02 

13.000000 

0.272152E-02 

14.000000 

0.333604E-C2 

15.0o0000 

0.398252E-02 

16.000000 

0.465989E-02 

17.000000 

0.536785E-02 

18.000000 

0.610650E-02 

19.000000 

0.687617E-02 

20.000000 

0.767725E-02 

21.000000 

0.851013E-02 

22.000000 

0.937520E-02 

23.000000 

0.102728E-01 

24.000000 

0.112031E-01 

25.000000 

0.121664E-01 

26.000000 

0.131629E-01 

27.0oOOOO 

0.141926E-01 

28.000000 

0.152556E-01 

29.000000 

0.16352OE-C1 

30.000000 

0.174816E-01 

31.000000 

0.186446E-01 

32.000000 

0.198408E-01 

33.0O0C 

0.21O70OE-01 

34.000000 

0.223322E-01 

35.000000 

0.236272E-01 

36.000000 

0.249549E-01 

37.000000 

0.263150E-01 

38.0o0000 

0.277073E-01 

39.000000 

0.291316E-01 

40.000000 

0.305878E-01 

41.000000 

0.320754E-01 

42*000000 

0.335943E-01 

43.0ououC 

0.351442E-01 

44.000000 

0.367248E-01 

45.000000 

0.383359E-01 

46.000000 

0.399771E-01 

47.000000 

0.416481E-01 

48.000000 

0.433486E-01 

49.000000 

0.450782E-01 

50.000000 

0.468368E-C1 

51.000000 

0.486239E-01 

52*000000 

0.504391E-01 

53.0ol,oOO 

0.522822E-01 

54.0ouuoo 

0.541528E-01 

55.000000 

0.560505E-01 

56.0o0000 

0.579750E-01 

57.000000 

0.599259E-01 

58.000000 

0.619029E-C1 

59.000000 

0.639055E-01 

60.000000 

0.659335E-C1 
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61.000000 

0.679864E-C1 

62.000000 

0.700638E-01 

63.000000 

0.721654E-01 

64.000000 

0.742908E-01 

65.000000 

0.764396E-01 

66.0OOOOC 

L.786114E-C1 

67.CuOOuO 

0.808059E-01 

68.000000 

0.830226E-01 

69.000000 

0.852610E-01 

70.000000 

0.875210E-01 

71.000000 

0.898020E-C 1 

72.000000 

0.921037E-01 

73.000000 

0.944256E-01 

74.000000 

0.967674E-C 1 

75.000000 

0.991286E-01 

76.000000 

0.101509 

77.000000 

0.103908 

78.000000 

0.106325 

79.000000 

0.108760 

80.000000 

0.111213 

81.000000 

0.113683 

82.000000 

0.116169 

83.000000 

C. 118672 

84.000000 

0.121190 

85.000000 

0.123725 

86.000000 

0.126274 

87.000000 

0.128838 

88.000000 

0.131416 

89.0UOOO0 

0.134008 

90.000000 

0.136614 

91.000000 

0.139233 

92.000000 

0.141865 

93.0U0OO0 

0.144509 

94.000000 

0.147166 

95.000000 

C. 149834 

96.000000 

0. 152513 

97.000000 

j. 155203 

98.000000 

0.157904 

99.000000 

0.160615 

100.000000 

0.163336 

101  .CouuUi, 

0.166066 

102.000000 

0.168806 

103.0UOO00 

0.171554 

104.00000C 

0.174310 

105.000000 

0.177075 

106.000000 

0.179848 

HI.  000000 

0.182627 

108.000000 

0.185414 

109.000000 

0.188208 

110.000000 

0.191008 

111.000000 

0.193814 

112.000000 

0.196625 
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113.000000 

0.199443 

114.000000 

0.202265 

115.000000 

0.205092 

116.000000 

0.207923 

117.000000 

0.210759 

118.000000 

0.213598 

119.0000PC 

0.216441 

120.000000 

0.219287 

121.000000 

0.222136 

122.000000 

0.224988 

123.0^0000 

0.227842 

124.000000 

0.230698 

125.000000 

0.233557 

126.000000 

0.236416 

127.000000 

0.239277 

128.000000 

0.242139 

129.000000 

0.245001 

130.000000 

0.247865 

131.000000 

0.250728 

132.000000 

0.253591 

133.000000 

0.256454 

134.000000 

0.259317 

135.000000 

0.262178 

136.000000 

0.265039 

137.000000 

0.267898 

138.000000 

0.270756 

139.000000 

0.273613 

140.00000C 

0.276467 

141.000000 

0.279319 

142.000000 

0.282169 

143.000000 

0.285017 

144.000000 

0.287861 

145.000000 

0.290703 

146.000000 

0.293542 

147.000000 

0.296377 

148.000000 

0.299208 

149.000000 

0.302036 

150.0OU000 

0.304860 

151.000000 

0.307680 

152.000000 

0.310496 

153.000000 

0.313307 

154.000000 

0.316113 

155.000000 

0.318915 

156.000000 

0.321712 

157.000000 

0.324503 

158.0O0000 

0.327290 

159.0O00U0 

0.330071 

160.0u0000 

0.332846 

161.000000 

0.335616 

162.0^0000 

0.338379 

163.000000 

0.341137 

164.000000 

0.343889 
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165.000000 

0.346634 

166.000000 

0.349373 

167.000000 

0.352106 

168.000000 

0.354832 

169.000000 

0.357551 

170.0o00o0 

0.360263 

171.0UOOOO 

0.362968 

172.0OO000 

0.365666 

173.000000 

0.368357 

174.000000 

0.371040 

175.000000 

0.373716 

176.000000 

0.376385 

177.000000 

0.379046 

178.000000 

0.381699 

179.000000 

0.384344 

180.000000 

0.386982 

181.000000 

0.389611 

182.0o0000 

0.392232 

183.0uwU00 

0.394846 

184.0QO000 

0.397451 

185.0000CO 

0.400047 

186.000000 

0.402635 

187.0OC00O 

0.405215 

188.000000 

0.407786 

189.000000 

0.410  349 

190.000000 

0.412903 

191.000000 

0.415448 

192.000000 

0.417985 

193.0OUO00 

0.420512 

194.000000 

0.423031 

195.000000 

0.425541 

196.000000 

C. 428042 

197.000000 

0.430533 

198.000000 

0.433016 

199.000000 

0.435489 

200.000000 

0.437954 

201.000000 

0.440409 

2C2. 000000 

0.442854 

203.000000 

0.445291 

204.000000 

C. 447718 

205.0O0000 

0.450136 

206.0O0000 

0. 452544 

2C7.0o000o 

0.454943 

208.000000 

0.457332 

209.000000 

0.459712 

210.000000 

0.462083 

211.000000 

0.464443 

212.000000 

0.466795 

213.000000 

0.469137 

214.000000 

0.471469 

215.0O0000 

0.473791 

216.000000 

0.476104 
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217.000000 

0.478407 

218.OUO0O0 

0.480701 

219.OU0OOO 

0.482985 

220.000000 

0.485259 

221.000000 

0.487523 

222.000000 

C. 489778 

223.000000 

0.492023 

224.000000 

0.494259 

225.000000 

0.496484 

226.000000 

0.498700 

227.000000 

0.500906 

228.000000 

0.503103 

229.0o0uuo 

0.505290 

230.000000 

0.507467 

231.000000 

0.509634 

232.000000 

0.511792 

233.000000 

0.513940 

234,000080 

0.516078 

235.000000 

0.518206 

236.000000 

0.520325 

237.000000 

0.522434 

238.000008 

0.524534 

239.0000OU 

0.526624 

24o.0uoouu 

0.528704 

241.000000 

0.530774 

2A2.quoo.OQ 

0.532835 

243.000000 

0.534887 

244.000000 

0.536928 

245.000000 

0.538961 

246.000000 

0.54C983 

247.000000 

0.542996 

248.000000 

0.545000 

249.000000 

0.546994 

250.000000 

0.548979 

251.000000 

0.550954 

252.0UOOUO 

0.552919 

253.000000 

0.554876 

254.0uO000 

0.556822 

255.0^0000 

0.558760 

256.0oouu0 

0.560688 

257.000ouO 

0.562607 

258.000000 

0.564516 

259.000000 

0.566417 

260.000000 

0.568308 

261.000000 

0.570189 

262.000000 

0.572062 

263.000000 

0*573925 

264.000000 

0.575779 

265.0^0000 

0.577624 

266.0^0000 

0.579460 

267.0u0uo0 

0.581287 

268.000000 

0.583105 

269.000000 

0.584914 
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270.000000 

0.586714 

771.000000 

0.588505 

772.000000 

0.590287 

273.000000 

0.592060 

274.000000 

0.593824 

275.000000 

0.595580 

276.000000 

0.597326 

277.000000 

0.599064 

278.000000 

0.600793 

279.000 

0.602514 

280.000000 

0.604226 

281.000000 

0.605929 

282.000000 

0.607623 

283.000000 

0.609309 

284.000000 

0.61C987 

285.000000 

0.612656 

286.0u0000 

0.614316 

287.000000 

0.615968 

288.000000 

0.617612 

289.000000 

0.619247 

290.000000 

0.620874 

291.000000 

0.622493 

292.0u0000 

0.624103 

293.000000 

0.625705 

294.000000 

0.627299 

295.000000 

0.628885 

296.000000 

0.630463 

297.000000 

0.632032 

298.0U0000 

0.633594 

299.000000 

0.635147 

300.000000 

0.636693 

301.000000 

0.638231 

302.0UOOOO 

0.639760 

303.000000 

0.641282 

304.000000 

0.642796 

305.000000 

0.644302 

306.000000 

0.645801 

307.000000 

0.647291 

308.000000 

0.648774 

309.000000 

0.650250 

310.000000 

0.651718 

311.000000 

0.653178 

312.0U0000 

0.654630 

313.0000^0 

0.656075 

314.000000 

0.657513 

315.000000 

0.658943 

316.0^0000 

0.660366 

317.000000 

0.661781 

318.000000 

0.663189 

319.000000 

0.664590 

320.000000 

C. 665984 

321.00000C 

0.667370 
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322.000000 

0.668749 

323.000000 

0.670121 

324.0^0000 

0.671486 

325*00001 

0.672844 

326.00C000 

0.674194 

327.000000 

0.675538 

328.000000 

0.676875 

329.000000 

0.678205 

330.000000 

0.679527 

331.000000 

0.680844 

332.000000 

C. 682153 

333.000000 

0.683455 

334.000000 

0.684751 

335.000000 

0.686040 

336.000000 

0.687322 

337.000000 

0.688597 

338.000000 

0.689866 

339.000000 

0.691129 

340.000000 

0.692384 

341.000000 

0.693634 

342.000000 

0.694876 

343.000000 

0.696113 

344.000000 

0.697343 

345.000000 

0.698566 

3*6.000000 

0.699784 

347.0u00u0 

0.700995 

348.000000 

0.702199 

349.000000 

0.703397 

350.000000 

0.704590 

351.000000 

0.705776 

352.000000 

0.706956 

353.000000 

0.708129 

354.000000 

0.709297 

355.000000 

0.710459 

356.000000 

0.711614 

357.000000 

0.712764 

358.000000 

0.713908 

359.000000 

0.715045 

360.000000 

0.716177 

361.000000 

0.717303 

362.000000 

0.718424 

363.000000 

0.719538 

364.000000 

0.720647 

365.000000 

0.721750 

366.000000 

0.722847 

367.000000 

0.723939 

368.000000 

0.725025 

369.0O000C 

0.726106 

37U.0O0000 

0.727181 

371.000000 

0.728250 

372.000000 

0.729314 

373.000000 

0.730372 
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374.000000 

0.731425 

375.000000 

0.732473 

376.000000 

0.733515 

377.000000 

0.734552 

378.000000 

0.735583 

379.0UOOOO 

0.736610 

380.000000 

0.737630 

381.000000 

0.738646 

382.000000 

0.739657 

3tQ3>0^00'00 

0.740662 

384.000000 

0.741662 

385.000000 

0.742658 

386.000000 

0.743648 

387.000000 

0.744633 

388.000000 

0.745613 

389.000000 

0.746588 

390.000000 

0.747558 

391.000000 

0.748523 

392.0o00uu 

0.749483 

393.0uO0OC 

0.750438 

394.000000 

0.751389 

395.000000 

0.752335 

396.000000 

0.753275 

397.000000 

0.754211 

398.000000 

0.755143 

399.000000 

0.756069 

400. 000000 

0.756991 

401.000000 

0.757909 

402.000000 

0.758821 

403.000000 

0.759729 

404.OO00U0 

0.760633 

405.000000 

0.761532 

406.000000 

0.762426 

407.000000 

0.763316 

408.000000 

0.764201 

409.000000 

0.765082 

410.000000 

0.765959 

411.000000 

0.766831 

412.000000 

0.767699 

413.000000 

0.768562 

414.0O000C 

0.769421 

415.0O000C 

0.770276 

416.000000 

0.771126 

417.000000 

0.771973 

418.000000 

0.772815 

419.000000 

0.773652 

420.000000 

0.774486 

421.000000 

0.775315 

422.000000 

0.776141 

423.000000 

0.776962 

424.000000 

0.777779 

425.000000 

C. 778592 
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426.Uv/Ooou 

0.779402 

427.0ooOoO 

0.780207 

428.000000 

0.781008 

429.000000 

0.781805 

430.0o0000 

0.782598 

431.000000 

G. 783387 

432.000000 

0.784173 

433.000000 

0.784954 

434.000000 

0.785732 

435.000000 

0.786505 

436.000000 

0.787275 

437.0O0000 

0.788042 

438.0gougg 

0.788804 

439.000000 

0.789563 

440.000000 

0.790318 

441.000oOC 

0.791069 

442.000000 

0.791817 

443.0o0000 

0.792561 

444.000000 

0.793301 

445.000000 

0.794038 

446.000000 

0.794771 

"*47.0oU0OO 

0.795500 

448.0000  JO 

0.796226 

449.0O0000 

0.796949 

450.000000 

0.797668 

451.0oOOOO 

0.798383 

452.000000 

0.799095 

453.000000 

0.799804 

454.000000 

0.800509 

455.000000 

0.801211 

456.000000 

0.801909 

457.000000 

0.802604 

458.000000 

0.803295 

459.000000 

0.803984 

460.000000 

0.804668 

4*1.000000 

0.805350 

462.000000 

0.806029 

463.000000 

0.806704 

464.000000 

0.807375 

465.0ou0oo 

0.808044 

466.000000 

0.808710 

467.000000 

0.809372 

468.000000 

0.810031 

469.000000 

0.810687 

470.000000 

0.811340 

471.000000 

0.811989 

472.00000C 

0.812636 

473.000000 

0.813280 

474.0jOUoo 

0.813920 

475.0O00O0 

0.814557 

476.0o0000 

0.815192 

477.000000 

0.815823 
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478.000000  0.816452 

479.000000  0.817077 

48u.0u0000  0. 81770o 

481.0O0000  0.818319 

482.000000  0.818936 

483.000000  0.819550 

484.000000  0.820160 

485.000000  0.820768 

486.000000  0.821374 

487.000000  0.821976 

488.000000  0.822575 

489.000000  0.823172 

490.000000  0.823766 

491.0u00U0  0.824357 

492.000000  0.824945 

493.0vJ0O0o  0.825531 

494.000000  0.826114 

495.000000  0.826694 

496.000000  0.827271 

497.000000  0.827846 

498.000000  0.828418 

499.000000  0.828988 

500.000000  0.829554 
501.000000    •    0.830119 
STOP   END  OF  PROGRAM  AT   STATEMENT  00  6  +  CO  LINES 
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500.000000 

0.829554 

600.000000 

0.875128 

700.000000 

0.905112 

800.000000 

0.925698 

900.000000 

0.940361 

1000.000000 

0.951139 

1100.000000 

0.959273 

1200.000000 

0.965554 

1300.000000 

0.970499 

1400.000000 

0.974459 

1500.000000 

0.977678 

1600.000000 

0.980328 

1700.000000 

0.982535 

18C0.0^0u00 

0.984392 

190U.OUOO00 

0.985969 

2000.000000 

0.987320 

2100.000000 

0.988485 

2200.000000 

0.989497 

2300.000000 

0.990382 

2400.000000 

0.991160 

2500.000000 

0.991847 

2600.000000 

0.992458 

2700.000000 

0.993002 

28C0. 000000 

0.993490 

2900.000000 

0.993929 

3060.0600OC 

0.994324 

3100.000000 

0.994683 

3200.000000 

0.995008 

33uO.Oui 

0.995305 

3400.000000 

0.995576 

3500.000000 

0.995824 

3600.000000 

0.996052 

3700.000000 

0.996261 

3800.000000 

0.996455 

3900.000000 

0. 996634 

4000.000000 

0.996799 

4100.000000 

0.996953 

4200.000000 

0.997096 

4300.000000 

0.997229 

4400.000000 

0.997354 

4500.0000UO 

0.997470 

4600.000000 

0.997578 

4700.000000 

0.997680 

4800.000000 

0.997775 

4900.000000 

0.997865 

5000.000000 

0.997949 

5100.000000 

0.998029 

STOP   END  OF 

PROGRAM  AT   STATEMENT  00 

6+00  LINES 
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5  OU 0  . 0  ^  0  0  o  0 

0.997949 

6OOO.O0OOOO 

0.998575 

7000.000000 

0.998953 

80C0. 000000 

0.999198 

9000.000000 

0.999366 

10000.  Oi-. 

0.999486 

11000. 000000 

0.999576 

12CO0. 000000 

0.999643 

13000;000000 

0.999696 

14000. OCOOOO 

0.999738 

15000.000000 

0.999772 

160C0i00i 

0.999799 

17000.000000' 

0.999822 

18CC0. 000000 

0.999841 

19000. OuOOOO 

0.999858 

20000.000000 

0.999872 

21000.000000 

0.999884 

22000.000000 

0.999894 

23000.000000 

0.999903 

24000.000000 

0.99Q91 1 

25000.000000 

0.999918 

26000.000000 

0.999924 

27000.000000 

0.999930 

28000.00 

0.999934 

29000.00 

0.999939 

30000*000000 

0.999943 

31000. OOOOOO 

0.999947 

32000.000000 

0.999950 

33000.000000 

0.999953 

34000.001 

0.999956 

35000.000000 

0.999958 

36000.000000 

0.999960 

37000.000000 

0.999962 

38000.000000 

0.999964 

39000.0- 

0.999966 

40000.000000 

0.999968 

41000.000000 

0.999969 

42000.001 

0.999971 

43000.000000 

C. 999972 

44000.000000 

0.999973 

45000.000000 

0.999975 

46000.000000 

0.999976 

47000.000000 

0.999977 

48000.000000 

0.999978 

49000.000000 

0.999979 

50000.000000 

0.999979 

51000.000000 

0.999980 

52000.000000 

0.999981 

53000. OuOOOO 

0.999982 

54000. OoOOOO 

0.999982 

55000. 0^ 

0.999983 

56000.000000 

0.999984 
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RC  notch  filters  with  equal  amplitudes  and  unequal  ampli- 
tudes at  zero  and  infinite  frequencies  in  the  frequency  response 
curve  are  described  in  the  literature.   Their  main  applications 
are  in  feedback  amplifier  problems. 

Three  notch  filters  were  investigated  in  this  paper:   the 
Goldman  type  RC  notch  filter,  a  more  general  notch  filter  of 
the  Goldman  type,  and  a  still  more  complex  RC  notch  filter. 
The  parameter  values  for  the  existence  of  a  notch  frequency  were 
found  for  each  network  by  requiring  that  the  capacitance  have  a 
maximum  value  in  each  filter.   In  the  third  case  the  capacitance 
was  a  function  of  two  variables,  resulting  in  an  interesting 
optimization  problem. 

The  first  and  third  networks  were  found  to  be  equal  ampli- 
tudes filters,  whereas  the  second  was  found  to  be  an  unequal 
amplitudes  filter.   The  quality  factor,  Q,  for  the  third  net- 
work was  found  to  be  lower  than  that  of  the  Goldman  type  filter; 
thus  the  more  complex  network  features  a  wider  bandwidth  than 
the  simple  Goldman  type  notch  network. 


